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FOREWORD 


Some  results  of  this  report  have  been  briefly  presented  before  as 
a  lecture,  entitled,  "Theoretical  Calculation  of  the  Diffuser  Efficiency  of 
Supersonic  Wind  Tunnels  with  Free  Jet  Test  Section’*,  at  the  Heat  Transfer 
and  rluid  Mechanics  Institute,  Berkeley,  California,  in  June  igl+9  by  the 
author.  Dr.  Rudolf  Hermann.  (See  Reference  25) 

This  report  represents  a  much  more  detailed  analytical  study  and 
a  comprehensive  comparison  with  available  experiments.  Both  the  above 
mentioned  lecture  and  the  present  report  deal  with  the  special  case  of- 
pressure  equilibrium,  when  the  pressure  surrounding  the  free  jet  is  equal 
to  the  pressure  in  the  Laval  nozzle  exit  cross-section* 

Another  report  by  the  same  author,  not  yet  published,  deals  with 
the  general  case  of  pressure  non-equilibrium  and  is  entitled,  “Diffuser 
Efficiency  of  Free  Jet  Supersonic  Wind  Tunnels  at  Variable  Test  Chamber 
Pressure."  (See  Reference  26) 

The  author  wishes  to  express  his  thanks  to  Dr.  R.  H.  Mills  for 
his  steady  interest  in  these  problems  and  his  stimulating  discussions;  to 
r#  long-time  associate j  for  his  careful  and  extensive  numerical 

evaluation  of  the  equations  and  laborious  preparation  of  the  graphs;  to 
Mr.  H.  U.  Eckert  for  critical  study  of  the  draft  and  his  valuable  comments 
to  Messrs.  J.  A.  Spanogle  and  M,  o.  Lawson  for  their  help  in  editing  the 
English  text;  to  Miss  M.  Le  Comte  and  Miss  V,  Bowman  for  preparing  the 
typewritten  drafts  and  final  copy. 
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ABSTRACT 


In  the  wind  tunnel  arrangement  under  consideration,  the  air 
leaves  the  Laval  nozzle  as  a  free  jet  and  is  recaptured  by  the  diffuser, 
which  is  of  the  convergent-divergent  design.  A  theoretical  analysis  of 
the  flow  process  through  this  type  of  supersonic  wind  tunnel  is  presented 
and  the  diffuser  efficiency  is  calculated  for  the  case  of  equilibrium 
between  test  chamber  pressure  and  pressure  in  the  nozzle  exit,  assuming 
one-dimensional,  inviscous,  steady  flow. 

Using  the  basic  equations  of  continuity,  energy  and  momentum 
flux  through  a  bounding  surface,  an  exact  solution  of  the  problem  is 
obtained,  which  is  applicable  up  to  Mach  number  infinite. 

One  of  the  basic  results  is,  that  in  the  recapturing  zone  of 
the  diffuser  a  transition  occurs  from  supersonic  to  subsonic  flow,  which 
is  followed  by  an  acceleration  in  the  convergent  portion  up  to  sonic 
velocity  at  the  second  throat.  The  transition  is  not  a  normal  shock  and 
involves  a  total  pressure  loss  greater  than  that  of  a  normal  shock  at  the 
test  section  Mach  number,  A  Mathematical  solution  with  supersonic  velocity 
after  the  transition  process  has  no  physical  existence. 

A  comprehensive  comparison  of  the  analytical  results  with 
available  experiments  in  supersonic  wind  tunnels  up  to  Mach  numbe  •  I+.I4. 
regarding  diffuser  efficiency  and  second  throat  area  shows  good  agreement. 
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DIFFUSES  EFFICIENCT  AND  FLOW  PROCESS  OF  SUPERSONIC 
WOOD  TUNNELS  WITH  FREE  JET  TEST  SECTION 


Section  1  -  Introduction  -  Diffusers,  Flow  Process  and  Efficiency  Definitions, 
A.  Diffusers  in  Wind  Tunnels  with  Free  Jet  and  Closed  Test  Section 


Supersonic  wind  tunnels,  like  subsonic  wind  tunnels,  may  be  placed 
in  one  of  two  categories  according  to  the  bounding  surface  of  the  air  stream 
in  the  test  section,  i.e.,  those  With  open  or  free  jet  test  section  and  those 
with  closed  or  rigid  wall  test  section.  The  characteristic  of  the  first  type 
is  the  constant  pressure  around  and  along  the  surface  of  the  free  jet,  while 
the  direction  of  the  streamlines  at  the  surface  is  not  determined.  The 
characteristic  of  the  second  type  is  the  fixed  direction  of  the  velocity  along 
the  rigid  walls  of  the  test  section,  while  the  pressure  inside  the  flow  or 
along  the  walls  of  the  test  section  are  dependent  variables.  It  is  well  known 
for  the  subsonic  incompressible  and  compressible  velocity  range,  that  the  in¬ 
fluence  of  the  limited  size  of  a  wind  tunnel  stream  on  the  aerodynamic  test 
data  (drag,  lift)  depends  on  the  boundary  conditions  and  the  so-called  wind 
tunnel  corrections  are  quite  different  for  the  two  types. 

A  similar  behavior  exists  with  regard  to  the  diffuser  efficiency  of 
a  supersonic  wind  tunnel.  The  following  analysis,  assuming  one  dimensional 
flow,  will  show  that  the  flow  process  in  a  diffuser  of  a  supersonic  wind 
tunnel  with  free  jet  test  section  is  completely  determined  by  the  preceding 
free  jet  and  its  spreading  due  to  the  turbulent  mixing  zones  at  its  surface. 

Thus  the  area  F2  at  the  entrance  of  the  supersonic  diffuser  must  be  larger 
than  the  cross  section  area  F^  of  the  stream  leaving  the  nozzle  (F2  >  F-, ). 

The  cross  section  ratio  f  -  Ft/Fo  enters  the  fundamental  equations  (Section  2), 
and  determines  its  solution  (Section  3).  It  is  therefore  an  important  para¬ 
meter  of  the  diffuser  efficiency  (Section  i+)  and  other  related  wind  tunnel 
characteristics,  for  instance,  the  second  throat  area  (Section  5)»  In  case 
of  a  supersonic  wind  tunnel  with  closed  test  section  this  parameter  does  not, 
of  course,  appear,  and  the  flow  process  in  the  diffuser  will  be  quite  different. 
Therefore,  it  is  necessary  to  differentiate  clearly  between  the  diffuser  effi¬ 
ciency  of  the  supersonic  wind  tunnel  with  free  jet  test  section  and  that  for 
the  closed  test  section. 

A  third  type  of  wind  tunnel  test  section  is  sometimes  used  or  proposed, 
the  so-called  ”half -open,  half -closed”  wind  tunnel,  which  offers  advantages 
in  minimizing  the  wind  tunnel  corrections  in  the  incompressible  and  in  the  com¬ 
pressible  subsonic  region.  Assuming  a  rectangular  shaped  cross  section,  two  of 
the  air  stream  surfaces  are  bounded  by  still  air  as  in  the  free  jet  type,  while 
the  other  two  are  bounded  by  rigid  walls  as  in  the  closed  type.  Analytical 
explanation  of  the  flow  process  in  the  diffuser  downstream  of  such  a  test  section 
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seems  to  be  more  difficult  than  in  the  case  of  the  two  other  types. 

It  might  be  briefly  mentioned  that  the  results  obtained  for  the 
diffuser  of  a  free  jet  supersonic  wind  tunnel  are  not  applicable  to 
diffusers  in  supersonic  ram  jets,  or  vice  versa.  The  basic  differences 
in  these  cases  are  as  follows: 

1.  Fixed  mass  flow,  given  through  the  preceding  Laval  nozzle,  in 
the  wind  tunnel  and  variable  mass  flow  with  traveling  shock  ("spilling") 
in  case  of  ram  jet. 

2.  Finite  size  of  the  jet  with  a  spreading  turbulent  mixing  zone 
collected  by  the  wind  tunnel  diffuser  and  infinite  uniform  medium  with¬ 
out  any  friction  layer  approaching  the  intake  of  the  ram  jet  diffuser. 

B.  Flow  Process  and  Tasks  of  the  Diffuser 


The  following  analysis  will  consider  a  conventional  type  of  a 
supersonic  wind  tunnel  with  a  free  jet  test  section  (See  Fig.  1)  -  to  the 
author’s  knowledge  first  established  in  Goettingen  by  L.  Prandtl  and  A. 
Busemann  in  1$J22>  (Ref.  1),  The  air  flows  from  the  open  atmosphere 
or  from  a  settling  chamber  through  a  straightener  and  entrance  cone  into 
the  Laval  nozzle  and  attains  sonic  velocity  in  the  minimum  section  called 
the  first  throat.  The  stream  leaves  the  exit  cross  section  of  the  Laval 
nozzle  with  supersonic  velocity  as  a  free  jet  bounded  by  air.  In  most  of 
the  supersonic  tunnels  an  air  tight  test  chamber  surrounds  the  free  jet, 
whose  pressure  then  is  independent  of  the  atmosphere.  The  free  jet  is 
collected  by  the  diffuser*  In  contrast  to  subsonic  wind  tunnels,  the 
supersonic  wind  tunnel  diffuser  consists  of  two  parts j  the  first  is  a 
converging  diffuser  portion  which  ends  up  in  a  minimum  section  called  the 
second  throat,  and  the  second  a  diverging  diffuser  portion,  which  finally 
leads  to  the  blowers. 

The  main  tales  of  the  diffuser  is  the  transfer  of  the  large  kinetic 
energy  at  low  pressure  of  the  test  section  stream  into  h igh  pressure  at 
nearly  vanishing  kinetic  energy  at  the  end  of  the  diffuser.  This  purpose 
is  the  same  in  all  subsonic  and  supersonic  wind  tunnels  with  either  free 
jet  or  closed  test  section.  Due  to  the  increasing  kinetic  energy  in 
supersonic  tunnels  and  ths  large  driving  power  needed  for  operation,  the 
efficiency  of  the  pressure  recovery  is  becoming  more  and  more  important. 

The  second  task  of  the  diffuser,  which  is  confined  only  to  supersonic  wind 
tunnels  with  free  jet  test  section,  is  the  control  of  the  air  pressure  pc 
in  the  test  chamber,  which,  at  supersonic  velocity  of  the  jet,  is  not  ^ 
necessarily  equal  to  the  pressure  of  the  air  stream  in  the  exit  cross 
section  of  the  Laval  nozzle.  To  fulfill  these  two  requirements,  the 
diffuser  is  usually  designed  with  adjustable  throat  cross  section  and 
adjustable  walls,  upstream  and  downstream  from  the  throat,  forming  a 
reasonable  aerodynamic  shape. 
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The  reason  for  building  the  first  part  of  the  diffuser  convergent  is  the 
assumption  more  or  less  generally  made  up  to  the  present  time  that  the  in¬ 
coming  supersonic  velocity  should  be  continuously  decelerated  which  results 
in  increasing  pressure.  The  optimum  recovery  would  be  reached,  if  the  sonic 
velocity  is  attained  in  the  second  throat,  so  that  subsonic  velocities  only, 
without  further  shocks,  occur  in  the  divergent  part  of  the  diffuser#  Actually, 
this  ideal  case  is  only  approached  as  the  velocity  in  the  throat  is  still 
slightly  supersonic  and  is  followed  try  a  normal  shock  close  downstream,  which 
leads  to  subsonic  velocity.  In  any  case,  the  convergent  part  of  the  diffusa: 
should  act  as  a  ■converted"  Laval  nozzle  declerating  the  supersonic  flow 
of  the  diffuser  entrance  to  approximately  sonic  velocity  in  the  throat  and 
to  subsonic  in  the  diverging  part0  It  follows  from  the  theoretical  one- 
dimensional  calculations  below  that  this  concept  is  wrong  for  the  open  jet 
tunnel#  When  the  jet  is  recaptured  a  sudden  transition  from  the  supersonic 
to  subsonic  flow  occurs  dose  to  the  entrance  section  of  the  diffuser  (See 
Sec.  I4C,  2  and  3),  which  is  followed  by  a  new  acceleration  in  the  converging 
part  of  the  diffuser  to  exact  sonic  speed  in  the  second  throat.  The  so-called 
supersonic  diffuser  acts  as  a  "conventional"  accelerating  Laval  nozzle,  not 
as  a  decelerating  "converted"  Laval  nozzle*  (See  Sec.  5A). 

The  mechanism  of  controlling  the  test  chamber  pressure  by  the  adjustable 
diffuser  throat  has  never  been  exactly  understood.  Actually  it  works  satis¬ 
factorily  in  many  but  not  in  all  cases.  Because  this  pressure  control  is  an 
"upstream  action"  in  supersonic  flow,  sometimes  the  thick  boundary  layers 
inside  the  diffuser  were  assumed  to  transmit  the  pressure  upstream.  This 
concept  is  incorrect  and  has  to  be  revised.  Pressure  control,  so  important 
for  actual  model  testing  in  the  tunnel,  can  be  explained  by  one-dimensional 
calculation  without  consideration  of  any  boundary  layer  (see  Sec.  5^). 

The  subsequent  calculations  were  instigated  by  some  results  of  measure¬ 
ments  carried  out  in  the  Peenemuende  supersonic  free  jet  wind  tunnel  by  my 
former  associate,  H.  Ramm  (Ref.  2)  in  1943*  reported  1945*  Pressure  dis¬ 
tribution  measurements  along  the  entire  wind  tunnel  and  especially  those  in 
the  diffuser  proved  that  most  of  the  losses  occur  in  the  first  part  of  the 
diffuser  and  are  not  distributed  over  the  whole  diffuser  and  tunnel.  There¬ 
fore,  it  seemed  logical  that  a  theoretical  analysis  of  the  flow  process 
between  nozzle  exit  and  diffuser  throat  would  produce  significant  results. 
Having  obtained  the  main  theoretical  results  in  1946  some  of  the  previous 
tests  results  acquired  greater  significance.  Afterwards  (in  1947) »  a  new 
detailed  discussion  of  his  measurements  was  given  by  H.  Ramm  (Ref.  3)* 

Nevertheless,  the  following  analysis  of  the  diffuser  efficiency  is 
purely  theoretical  and  does  not  use  any  empirical  values.  A  comparison 
of  the  results  with  all  experimental  data  available  to  the  author  is  given 
in  Section  6. 

C.  Various  Definitions  of  Diffuser  Efficiency  and  Energy  Ratio 

To  obtain  a  convenient  measure  for  the  pressure  recovery  of  the 
diffuser  various  definitions  of  different  kinds  of  efficiency  have  been 
introduced.  If  they  are  defined  correctly,  that  means  utilizing  thermo¬ 
dynamic  states,  they  can  be  transformed  one  into  another  as  a  function 
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of  Mach  No,  and  it  is  a  matter  of  convenience  which  one  to  use.  Some  of 
those  definitions  which  fulfill  the  mentioned  requirement  are  discussed 
in  the  following. 


Pressure  Efficiency  and  Energy  Efficienc 


J,  Ackeret  (Ref,  1|)  gives  two  possible  definitions  to  use  far 
evaluation  of  experimental  data  in  the  compressible  subsonic  region.  However, 
he  does  not  give  theoretical  calculations  of  any  diffuser  flow  process  or 
®ff iciency,  one  of  these  definitions  is  the  “pressure  efficiency1*  given  by 

V,  _  le  -  Pi 

P,  -  pi 

where  pQ  ■  total  pressure  before  the  diffuser,  pi  »  static  pressure  in  the 
inlet  cross  section  of  the  diffuser;  p<?  *  static  or  total  pressure  in 
the  exit  cross  section  of  the  diffuser  where  the  velocity  head  can  be 
neglected,  (For  evaluation  purposes  the  total  pressure  is  taken).  It 
compares  the  actual  pressure  increase  in  the  diffuser  with  the  maxi mum 
possible  increase  at  isentropic  flow.  Obviously  this  definition  originated 
in  the  theory  of  incompressible  flow,  where  the  denominator  represents  the 
kinetic  energy  in  the  smallest  section  of  a  venturi  meter.  The  pressure 
efficiency  for  reasonable  diffusers  in  the  incomp re ss ible  or  subsonic 
compressible  range  is  of  the  order  of  0,S,  At  high  supersonic  Mach  Nos,  it 
decreases  rapidly  and  drops  down  to  inconveniently  small  values;  far  instance, 
■  0,05  at  Mach  No,  5  and  ^p  “  0,002  at  Mach  No,  10, 


The  other  definition  is  the  " 


energy  efficiency"  given  by 

Jjl£!L 

To-7i 


where  T  *  stagnation  temperature;  Tt  “  temperature  in  the  inlet  cross  section 
of  the  aif fuser;  %  ■  blower  outlet  temperature  after  adiabatic  compression 
of  the  air  leaving  the  diffuser  up  to  the  original  total  pressure  p0«  It 
compares  the  difference  of  the  kinetic  energy  in  the  test  section  and  the 
blower  energy  with  t  he  kinetic  energy  in  the  test  section.  Obviously  this 
definition  is  desired  for  wind  tunnel  purposes,  where  a  complete  cycle  is 
required  and  the  blower  energy  necessary  to  maintain  this  cycle  is  of  pri¬ 
mary  interest.  The  “energy  efficiency*8  in  the  subsonic  range  for  reasonable 
diffusers  has  similar  values  (about  0»S)  as  the  "pressure  efficiency," 

Its  major  disadvantage,  however,  is  that  its  value  decreases  rapidly  with 
Mach  No,,  passes  zero  and  becomes  negative  at  higher  Mach  N0s.,  even  at 
diffusers  with  reasonable  pressure  recovery,  (Far  instance,  at  Ma  *  i+j 
**} t  “  “0,25;  at  Ma  ■  10:  "  -Au0)o  For  a  bad  diffuser  with  ho' pressure 

recovery  at  all,  the  energy  efficiency  becomes  zero  for  incompressible  flow 
only,  but  negative  even  at  subsonic  compressible  Mach  Nos,  (See  Sec,  5,  D)» 


2, 

L,  Crocco  (Ref,  5)  has  introduced  another  definition  of  the 
diffuser  efficiency,  assuming  a  step  by  step  compression  along  the  diffuser. 
He  defines  differentially  an  efficiency  T^at  the  ratio  between  actual  and 
isentropic  relative  increase  of  static  pressure  with  correspondingly  equal 
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decrease  of  kinetic  energy*  ^Cr  is  unknown 

the  whole  compression*  This  way  he  obtains 

tions  used  above,  can  be  transformed  into 
* 


^Cr  Jm. 


1L 


_e» 

Pi 


but  is  assumed  constant  over 
a  formula,  which  in  the  nota- 


Crocco  gives  no  theoretical  calculation  of  the  diffuser  process  or  the 
efficiency  itself.  He  calculates  with  the  aid  of  this  formula  by 
evaluation  of  total  pressure  measurements  pe  of  specific  wind  tunnel 
tests*  These  efficiency  values  converted  into  the  definition  used  in 
this  paper  are  discussed  later  (see  Sec.  6,  E  and  Fig*  20*  Curve  6)* 


Efficiency  values  ^ of  reasonable  diffusers  avoid  the  above  men¬ 
tioned  inconveniences  involved  in  Ackeret's  Vp  or  values  at  high  Mach  Nos0 
and  stay  in  a  convenient  order  of  magnitude  (for  instance,  0«lt  at  Mach  No* 

10)*  At  pressure  ratios  near  unity  the  formula  for  ^approaches  asympto¬ 
tically  the  formula  for  ,  thus  furnishing  conventional  values  of  the 
order  of  0.8  in  the  subsonic  range*  However,  Crocco  seems  not  to  differ¬ 
entiate  between  diffusers  behind  closed  or  free  jet  test  section »  In  case 
of  a  free  jet  tunnel,  his  concept  of  a  step  by  step  compression  along  the 
diffuser  is  incorrect  according  to  the  results  of  thB  following  analysis. 

Hence  the  values  of  efficiency  ^obtained  thereby  are  fictitious.  In  spite 
of  their  convenient  order  of  magnitude,  they  have  no  physical  meaning* 

3*  The  Energy  Ratio 

In  order  to  describe  the  performance  of  wind  tunnels  the  so- 
called  energy  ratio  is  also  often  used.  Hereby  the  kinetic  energy  inside 
the  test  section  1^.  is  considered  as  desired  energy  and  compared  with  the 
energy  consumption  of  the  biowers  Lg  LKin 

*  "*  Lb 

L  is  generally  defined  as  the  work  per  unit  time  either  at  isothermic  or 
at  adiabatic  compression,  and  £  is  called  the  isothermic  or  adiabatic  energy 
ratio  respectively*  Sometimes  the  energy  output  of  the  driving  motor  of  the 
blower  is  considered  as  L_*  A  detailed  discussion  of  the  values  of  the  energy 
ratio  calculated  with  the2  diffuser  efficiency  obtained  by  this  analysis  is 
given  later  (See  Sec*  5 j  B,  and  Fig,  16,  17) •  Tt  might  be  mentioned  that  the 
fraction  in  the  definition  of  the  energy  efficiency  is  the  reciprocal  value 
of  the  adiabatic  energy  ratio. 

In  discussing  the  definition  of  energy  ratio  we  have  to  consider  that 
the  kinetic  energy  in  the  test  section  1*^  is  desired  and  needed  to  fulfill 
the  task  of  the  tunnel.  But  a  wind  tunneiris  not  an  engine  which  takes  power 
L  from  outside  and  delivers  as  external  useful  work.  The  steady  flow 

o?  air  in  a  wind  tunnel  from  the  stagnation  chamber  to  the  section  upstream 
of  the  blower  has,  even  with  losses  in  the  diffuser  or  elsewhere,  constant 
energy,  if  all  kinds  of  energy  are  considered  and  non-conduting  tunnel  walls 
are  assumed.  In  case  of  wind  tunnels  with  a  closed  circuit  the  air  circulates 
in  a  complete  thermodynamic  cycle.  The  energy  Lg  put  into  the  air  flow  by 
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fch*  motor  of  tho  blower  ho*  to  bo  removed  entirely  a*  hoot  by  tho 
ooeling  ayttoa.  If  tho  compressor  pooor  i*  usod  In  defining  energy 
ratio  or  efficiency  of  the  diffu*er  or  of  the  whole  tunnel,  it  ie 
neeessary  to  presort be  a  specifio  compression  and  cooling  process. 

This  be  an  adiabatic  or  isothemio  or  any  other  kind  of  poly- 

tropio  compression  and  a  oooling  process  between  arbitrarily  defined  -values 
of  pressures  and  temperatures.  Cooling  can  take  place  before,  in  between, 
or  behind  the  compression.  Therefore,  the  definition  of  the  energy  ratio 
or  the  energy  efficiency  of  the  diffuser  depends  on  a  specifio  but 
arbitrarily  chosen  cycle. 


Efficiency  Definition  Psed  in  this  Paper 


The  definition  for  effioiency  used  in  this  paper  is  different 
from  those  mentioned  above.  All  so-called  losses,  which  occur  by  friction 
along  tha  tunnal  walls,  by  the  model  drag,  by  compression  shook*,  in  the 
diffuser  or  elsewhere  actually  mean  no  lose  but  a  devaluation  of  total 
energy,  vhich  ie  equivalent  to  an  increase  of  entropy.  From  this  point 
of  view,  the  air  flow  in  a  wind  tunnel  can  be  considered  in  the  same  manner 
as  any  compressible  flow  at  constant  energy  but  with  entropy  increase.  In 
order  to  characterize  such  an  entropy  increase  for  aerodynamio  purposes 
it  is  convenient  to  use  the  ratio  of  the  total  pressures  p£  /  pQ  after 
and  before  the  losses  occur,  as  done  in  case  of  compression  shook#  by 
A.  Busemann  (Kef.  6,  Par.  27,  page  1+37).  The  total  pressure  ratio  was 
widely  used  also  by  *.  Tollmien  (Ref.  7)  and  called  "Drosselfaktor"  K 
•. throttling  factor).  It  is  connected  with  the  entropy  increase  from 
»o  to  a  o  by  th«  equation  , 


The  use  of  total  pressure  ratio  is  therefore  equivalent  to  the  determination 
of  the  entropy  increase;  however,  no  energy  consideration  is  involved. 


Defining  the  diffuser  efficiency  we  compare  the  total  pressure 
ratio  of  the  diffuser  Kj>  with  the  total  pressure  ratio  of  a  normal  shock 
5  ws  At  th«  test  section  Mach  fcusnbep 

%  J£l  *  L^L\  •  (JL) 

Kns  \  P*  mfrustR  *  \  P»  /normal  shock 

ft#  call  this  coefficient  of  total  pressure  ratios  for  tho  purpose  of 
briefness  "diffuser  efficiency",  even  though  it  is  not  a  ratio  of  energies 
and  its  value  can  be  smaller  or  larger  than  one.  (For  subsonic  test  section 
velocities  the  denominator  is  taken  as  unity). 


The  order  of  magnitude  of  the  values  actually  obtained  will  show 

definition,  especially  in  the  range  of  high  Mach 
Numbers.  This  is  valid  for  wind  tunnel  diffusers  after  free  jet  or  closed 

:*!;JfCti0n  **  W®U  41  for  raa>t  diffusers.  Moreover,  this  efficiency 
definition  is  suggested  in  the  case  of  the  free  jet  wind  tunnel  diffuser 
by  the  results  of  the  following  analysis.  They  show  that  the  flow  process 
in  the  diffuser,  while  not  that  of  a  normal  shook  does  have  soma  similarity 
to  a  normal  shook.  Besides,  thsre  exists  a  finite  non-vanishing  limiting 
valua  of  %  for  Mach  So.  infinity,  which  is  approached  very  well  already 
ffom  Mach  Nos.  3  or  5  upwards.  None  of  the  other  definitions  mentioned 
above  offer  the  same  advantage.  This  definition  has  been  used  before  by 
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the  author  and  his  associates  at  Peenemuende  and  Kochel,  Germany,  since 
I9I42  for  wind  tunnel  design  purpose  without  knowledge  of  the  subsequent 
calculations. 


Section  2  -  Establishing  the  Fundamental  Equations 

A.  The  Spreading  of  the  Jet  and  the  Transformation  Zone  in  the  Diffuser 

The  airstream  emerges  from  the  exit  cross  section  Ft  of  the  Laval 
nozzle  (See  Fig.  1)  with  supersonic  velocity  w-,  as  a  free  jet  surrounded  by 
the  air  in  the  test  chamber  which  is  supposed  to  be  at  rest  in  a  greater 
distance  from  the  jet.  The  length  1  of  the  free  jet  depends  upon  various 
test  requirements  and  is  conventionally  between  approximately  one  and  two 
times  the  height  H  of  the  Laval  nozzle  exit.  The  free  jet  spreads  out,  in¬ 
creasing  its  area,  due  to  a  turbulent  mixing  process  at  its  surface,  resulting 
in  wedge  shaped  turbulent  mixing  zones.  This  process  has  been  analyzed  for 
incompressible  flow  by  W.  Tollmien  (Ref.  g)  and  A.  Kuethe  (Ref.  9)  using 
the  concept  of  the  turbulent  mixing  length  of  L.  Prandtl  and  recently  by 
H.  Liepmann  and  J«  Laufer  (Ref.  10)  with  detailed  results  on  the  turbulent 
structure  in  accordance  with  von  Karman's  (Ref.  11)  concept.  To  obtain  any 
pressure  recovery  the  free  jet  has  to  be  brought  into  the  diffuser.  This 
requires  a  "recapture"  of  the  free  jet  in  the  intake  Fg  of  the  diffuser.  Here 
the  free  jet,  with  a  given  pressure  as  the  boundary  condition,  is  transformed 
into  a  flow  between  rigid  walls  with  a  given  direction  of  velocity  as  determined 
by  the  walls.  We  assume  that  in  this  region  a  complete  transfer  of  all  values 
of  state  occurs  which  results  in  constant  velocity,  pressure  and  density  all 
over  the  cross  section.  Furthermore,  we  assume  that  this  transformation  takes 
place  ip  the  intake  portion  of  the  diffuser  with  constant  or  nearly  constant 
cross  sectional  area,  corresponding  to  a  slender  diffuser  design  as  used  in  many 
free  jet  supersonic  tunnels,  i.e. ,  in  Goettingen,  Aachen,  Peenemuende,  and 
Kochel  (See  Sec.  6,  D). 

The  spreading  of  the  free  jet  requires  that  the  intake  cross  section 
of  the  diffuser  F2  should  be  larger  than  the  Laval  nozzle  exit  cross  section  F^0 
This  is  a  conventional  arrangement  in  all  free  jet  wind  tunnel  designs.  (See 
Fig.  lg  and  Sec.  6,  A).  Often,  the  entrance  cross  section  of  the  diffuser  is 
equipped  with  a  bell-mouth,  a  knife  edge,  or  another  differently  shaped  pick-up. 
Such  a  device  is,  of  course,  important  to  the  separation  of  the  induced  air 
from  the  nozzle  air.  Its  interaction  with  the  mixing  zone  and  the  determination 
of  the  minimum  area  Fg >  which  is  required  for  a  smooth  recapturing  process  can 
be  analyzed  but  need  not  be  discussed  in  the  present  report.  Here  F2  is  one 
of  the  given  boundaries. 

B.  Test  Chamber  Pressure 

In  case  of  a  jet  with  supersonic  velocity,  the  test  chamber  pressure 
Pc  is  in  general  not  equal  to  the  pressure  p^_  in  the  nozzle  exit,  pcjt  p^.  Only 
a  specific  opening  of  the  second  throat  will  produce  pressure  equilibrium  pc  ■  p^. 
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corresponding  to  a  parallel  jet  in  the  core  inside  of  the  turbulent  mixing 
zones.  Larger  second  throat  areas  result  in  p  <  P-,,  smaller  ones  in  pQ  >  P, 

Test  chamber  pressure  pQ  >  Px  can  be  operated  but,Jlimited  by  the  requirement 
that  only  oblique  shocks  occur  at  the  nozzle  exit,  results  in  a  reduced  size 
of  test  rhombus  with  parallel  flow.  If  p  is  so  large  that  a  normal  shock 
occurs ,  no  supersonic  jet  can  be  maintained.  In  case  of  a  subsonic  jet 
pressure  equilibrium  pc  *  p]_  is  established  automatically,  and  the  second 
throat  area  determines  the  velocity  in  the  test  section. 

The  present  analysis  is  made  for  the  special  case  of  pressure  equi— 
librium.  A  solution  for  the  case  of  pressure  non— equilibrium  is  given  in 
another  report  (Ref.  26).  It  clarifies  the  exact  mechanism  of  test  chamber 
pressure  control  by  means  of  an  adjustable  diffuser  throat  and  gives  a  quan— 
titative  explanation  of  a  slew  starting  process  of  the  tunnel. 

Co  Assumptions  and  Fundamental  Equations 

We  consider  the  flow  as  one-dimensional  and  steady  and  use  the 
theorem  of  momentum  for  the  axial  momentum  flux  over  a  rectangular  shaped 
bounding  surface  which  extends  from  the  nozzle  exit  to  the  end  of  the  trans¬ 
formation  zone  and  has  a  constant  cross  section  Fg.  The  great  advantage  of 
such  application  of  the  momentum  theorem  is  that  the  internal  processes,  such 
as  mixing  at  the  jet  surface,  recapturing  and  transformation  in  the  diffuser 
intake  section,  need  not  be  known  in  detail.  The  final  state  w„,  p3,  9  P  in 
?2  after  completion  of  the  transformation  can  be  calculated  by  the  ^initial  state 

W>  5  In  and  the  amount  of  the  boundaries  F-,  and  Fp. 

Evidently  the  whole  transformation  process  is  irreversible  and  is  connected 
with  an  entropy  increase  which  will  eventually  be  determined  by  the  calcula¬ 
tion* 


The  following  detailed  assumptions  and  simplifications  are  made: 

1.  Steady  one-dimensional  flow  parallel  to  the  axis  is  assumed. 

This  simplification  is  more  valid  when  slender  diffusers  without  sudden 
variation  in  cross  section  are  considered.  Especially  the  intake  portion  of 
the  diffuser  must  be  long  enough  to  guarantee  a  complete  transformation 
process  (Sec,  2,  A).  ' 

2®  Air  tight  test  chamber  and  chamber  pressure  constant  with  time 
are  assumed.  The  air  induced  by  the  jet  out  of  the  test  chamber  must  be 
separated  and  returned  to  the  test  chamber,  before  the  transformation  process 
as  completed.  This  results  in  equality  of  mass  entering  and  leaving  the  bounded 
volume  and  thus,  the  continuity  equation  becomes 


(la) 


AFTR-633U 


8 


The  ••sumption  of  an  air  tight  t««t  chamber  is  of  greater  importance 
than  might  be  realised.  Calculations  show  that  relatively  small  in¬ 
leakage  of  air  into  the  test  section  oan  produoe  quite  a  different  flow 
prooess. 

3.  The  flu*  of  axial  momentum  through  the  bounding  surfaces  parallel 
to  the  axis  and  through  the  surface  normal  to  the  axis  in  the  upstream 
section  outside  Fi  is  negleoted,  because  the  components  of  the  induoed 
velocities  are  small  compared  with  w,.  Thus,  F.  ?,  w,1  is  the  momentum 
entering  and  Fz  w £  the  momentum  leaving  the  volume  within  the  bounding 
surface.  This  is  valid  for  a  test  chamber  of  reasonable  sis#  *iieh  allows 
the  turbulent  mixing  tones  to  develop  in  a  similar  manner  as  it  is  known  for  jets 
in  indefinite  medium.  Very  small  site  of  the  test  chamber  might  produce  addi¬ 
tional  Induoed  veleoities. 

I4.  The  statlf  pressure  in  the  test  chamber  pc  is  assumed  to  be  constant 
along  tha  bounding  surface  normal  to  the  axis  in  the  upstream  cross  saotion 
outside  F,,  which  is  identioal  with  neglecting  the  dynamic  pressure  of  the 
induoed  velocities.  (See  Par.  3)*  Henes,  the  foreei  aotlng  at  the  upstream 
cross  section  are  F,p,  +  (Fi- F, ) p<;  while  the  foroes  aoting 

at  tha  downstream  oross  saotion  era  F2  pz 

5*  The  friction  along  the  rigid  walls  Inside  the  diffuser  and  the 
turbulent  shearing  stresses  along  the  bounding  surfaoe  outside  the  diffuser 
are  neglected*  Therefore,  no  other  foroes  than  pressure  forces  according 
to  paragraph  U  enter  the  momentum  equation.  This  equation  states  that  the 
increase  in  momentum  of  the  gas  per  unit  time  equals  the  net  fore#  aoting 
on  the  gas  in  the  same  direction,  and  therefore  heeometi 

Fz  Wi*  —  Fl  9,  W,2  a  F,p,  (Fi-F,)pc  -  F2p2  (lb) 

It  ia  valid  for  the  general  case  of  pressure  non-equilibrium.  A  more  de¬ 
tailed  Investigation  concerning  the  effeet  of  the  turbulent  shearing  stresses 
along  the  jet  surfaoe  on  the  momentum  equation  (lb)  and  the  energy  equation 
(lo)  shows  that  tha  assumptions  made  in  Par.  5  and  6  hold  much  better  than 
might  bs  anticipated* 

6.  The  wind  tunnel  walla  are  assumed  to  be  Insulated,  so  that  no 
heat  is  transmitted  out  of  the  interior.  So  exterior  work  is  performed  by 
the  gas  avan  if  tha  work  of  the  turbulent  shearing  stresses  along  the  jet  surfaoe 
on  the  oh amber  air  is  taken  into  aoeount.  In  ease  of  an  airtight  test  chamber 
tha  air  within  belongs  to  the  system  because  it  is  in  permanent  energy  ex¬ 
change  with  the  stream  emerging  from  the  nos  tie.  Henoe  the  total  energy  (kinetic 
energy  plus  enthalpy)  entering  the  bounded  volume  must  equal  the  outgoing  total 
energy,  even  though  dissipation  is  present.  Thus,  the  energy  equation  gives 

_^L+CpT,  =-^-+CpT2«CpTo 

^vsrecpT^  Is  ths  enthalpy  of  the  air  at  rest  before  entering  the  tunnel. 
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Using  the  equation  of  state: 

this  is  transformed  into 

w,2  .  i  P, 

—  +  »-,  T, 

Thus,  we  ha-ve  a  set  of  3  independent  equations  (la,  lb,  lc)  for  the  solution 
of  the  three  unknown  variables  p,  9*  }  if  the“varxables  w,  p,  9, 
the  test  chamber  pressure  p  and  the  boundary  areas  Ft  and  are  given.  The 
system  of  equations  (la,  lb,  lc)  obtained  under  the  given  assumptions  for 
the  free  jet  diffuser  process  is  obviously  identical  with  that  of  a  com¬ 
pressible  pipe  flow  with  sudden  increase  in  cross  section  (See  for  instance  - 
A.  Busemann,  Ref.  6,  page  1*03).  However,  no  detailed  solution  and  discussion 
in  relation  to  the  free  jet  wind  tunnel  diffuser  problem  as  analyzed  here 
is— to  the  author’s  knowledge — given  elsewhere.  The  losses  connected  with 
the  pipe  flow  are  sometimes  called  "Carnot  impact  or  shock  loss”,  although 
impact  does  not  enter  into  this  problem  and  the  term  "shock”  has  no  relation 
to  the  compression  shock  waves  at  supersonic  flow,’ 


E  =T?T 

9 

-  .  1  Pa 

2  *-■  92  (lc) 


Section  3  -  Solution  of  the  System  of  Equations 

In  case  of  pressure  equilibrium  in  the  test  chamber  pc  *  pj_,  the  system 
of  the  continuity  equation  (la)  and  momentum  equation  (lb)  becomes: 

F.  vv,-  gt  *“f2-vv2.  92 

~  F,9(<  »  h(P,-p> 


(2) 


Combining  both  equations  gives 

F,w,  S/wj-w.)  =  F2(p,-p3) 

Using  the  abbreviation  ■—  *  f  <  !  one  obtains 

rz 

-E'—P*  »fW,.  9,  > 

Wi— J  I  JI  1 

-#1  9* _ f  f 

*>  §.  7  J 

W®  introduce  the  following  dimensionless  values  for  velocity,  pressure,  and 
mass  flow:*  '  9 

-X~=.  M  .  p  -  -  P  .  w  SL 

a*  '  >  po  -*■  5 

and  using 


(2a) 


0 


(2b) 


♦in* 


a*1? 


y~i 


^Footnote:  Distinguish  M  «  and  Ma 
denote  Mach  No. 


Mach  No.j  sometimes  M  is  used  to 
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«*  obtain  th*  relations 


©2  *  f-(9, 


System  (3)  plus  th*  energy  theorem  (ljs)  converted  into  th*  mum 
dimensionless  values  determines  the  three  unknown  variable*  Mz j  P*  j  ®i 

if  t  end  Mji  Px*  ©,  *r*  given. 

Th*  initial  *tat*  M,  Pi®,  ha*  th*  eta gnat ion  pressure  p0  and 
th*  initial  entropy  *j.  The  final  *tat*  Mg  Pg  ©^has  th* 

stagnation  pressure  p’  and  the  final  entropy  *g>  Me  introduce  the 
ratio  K  of  th*  total  pressures  before  and  after  th*  losses  ooour 

K  -  -£  (4a.) 

which  is  connected  with  th*  entropy  inorease  by  the  equation 


Sa-S. 


—In  K 


Then  w*  express  P  and  <9  as  funotlons  of  M  at  constant  entropy  s^ 
by  application  of  the  energy  equation  (le)  for  isentropio 
processes:  -  £ _  » 

®  (M;  s.)  -  f '  M  [  I  -  -5ir  M  •  J  J  (5) 

To  obtain  P  and  ©  as  funotions  of  M  for  another  constant  entropy 
Sg,  th*  right  hand  side  of  th*  equations  (5)  have  to  be  multiplied  by 
tne  total  pressure  ratio  K.  Or  in  other  words:  Pressure  and  mass  flow 
for  different  entropy  are,  at  any  velocity,  proportional  to  th* 
respective  stagnation  pressures.  For  this  reason  and  with  the  us*  of 
P  and  ®  as  function  symbols,  the  values  in  equation  (3)  may  be  written 

X  -  HM,;0  |  P,  =  K-T(M„s,)  \ 

®,  ■  |  ®,  =  -  K'®(M,i  S.)  J 


Equation  (3)  with  (5a)  becomes 

k-  P (yti±y)  -  2.  (MiiiA 

M,-Ma 


•f  •  0  (M,  iS.) 


K-@CM.is,)  -  )  J 

This  sytem  (6)  of  two  equations  determines  the  two  unknown  variables 
*2  *nd  K»  if  and  f  are  given. 
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Elimination  of  K  and  omission  of  tha  parameter  §,  finally  results 
in  one  equation  for  the  one  unknomn  Mpt 

*  r 

X OHil  _  21  _  V  f  1  \  £  (m  m  \ 

7  ®(K)  <S>(M.)  *W>77  '*'("»-"*] 

Equations  (5)  give  by  dividing  each  other 

-§-M  («) 

This  ratio  of  dimensionless  pressure  and  mass  flow  is  independent 
of  any  special  entropy  value.  The  transcendental  equation  (7)  la 
reduced  by  means  of  (g)  to  a  quadratic  equation  in  lig, 

f(7Tr)  "7^0  '  (jsr) 

This  is  simplified  finally  by  some  algebraic  transformations  toi 

7 -sr +(-Frr  -  t"£H  M' 


+-  M3 


or 


M‘‘- [f  T7  +  -  T  MJ Ml+ 1  -° 


(9) 


(9a) 


Tha  two  solutions  of  this  quadratic  equation  are  designated  by  M2ot  (+sign) 
and  Mip  (-signj 

( JA. - L  ilL.) M  %  i  -tjFf  it  l  t-t  W~  i  Ti  1 

Mg  -^r 1  f  <io) 


vefoci-Kes 


Hence,  the  final-«%&%&«  M  2&  can  be  calculated  by  (10),  if  the 
initial  swiffir  is  given*  p  It  follows  from  (9a)",  that 


M2*  •  M*ft  »  I 


(11) 


**r#  Mj.cc  >  I  and  Mifl  <  I 


Henoe,  it  should  be  emphasised,  that  our  problem  has  on*  supersonic 
solution  M2i  and  one  subsonio  solution  .  Their  velocities  are 

connected  in  the  same  way  as  tha  velocities  before  and  behind  the 
normal  shock.  It  is  demonstrated  later  (See  Sec.  >4,  A 2)  that  the 
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thermodynamio  states  (Mg 5 Pg)  also  of  both  solutions  srs  intsroonnsotsd 
by  a  normal  shook* 


It  must  bo  investigated  oarofully  whether  on#  or  both  solutions 
actually  axlst  in  aooordanoo  with  tho  seoond  lav  of  thormodynamioo  by 
determination  whether  tho  entropy  of  the  final  state  has  not  deoreased 
compared  with  that  of  the  initial  state  (See  Seo*  1»*  C)« 

of 

We  determine  now  the  second  unknown  equation  (6),  the  total  pressure 
ratio  I*  which  we  designate  from  now  on  as  Ip  because  it  denotes  the 
diffuser  losses.  Aeoording  to  the  two  solutions  for  M,  we  obtain  two 
values  of  Kp 


c  ©(M»i*.) 

T  ©(**)«.) 


^  K* 


r  ©  (  Mps,) 
T  ©(M^s.) 


(12) 


Now,  the  solution  of  the  problem  is  completed.  Equation  (10)  gives  Mg* 
(12)  gives  Ip  and  (5®)  with  (5)  determines  Pg  and  ©2 

Thus,  velocity,  pressure,  mass  flow  (and  therefore  density),  and  total 
pressure  ratio  of  the  final  state  are  know.  As  the  nest  step  the 
diffuser  efficiency  will  be  ealeulated. 


For  this  ealeulation  we  introduoe  by  definition  (See  Seo.  1,  C,  1*) 
the  diffuser  efficiency  as  the  total  pressure  ratio  of  the  diffueer  Ip 
compared  with  the  total  pressura  ratio  INS (K  i  )  of  a  normal  shook 
at  the  tost  section  velocity  Mj.  Ae  a  result,  it  follows 


yj  _ _ Kp«t  _  x  _ ©  (~Mi) _ __ 

t5flL_  K*is(M,)~  7  ©OW'Kns CM,) 

y,  =  Kpn  _  r  <9fM ,) 

K«s(Mi)  7  ©CM1(i)-  Knj(M,) 

Equations  (12)  and  (13)  allow  oaloulation  of  total  prasaura  ratio  and 
affioienoy  of  the  diffuser.  For  further  dlsoussion  and  numerical  cal¬ 
culation  of  ,  the  introduetion  of  a  new  funotlon 


®(M 

Kns(M) 


(1U) 


defined  in  the  supersonic  range  only  is  advantageous.  It  is  the  ratio 
of  the  narrowest  downstream  arose  seetlon,  attainable  after  a  normal 
shook  at  M  has  ocourred,  to  tho  shock  area.  KN$oan  he  written  ae  follows! 


K„SCM)  = 


1 


With  ©  aeoording  to  (5)>  ^  becomes 


w-tsT-K  i- 


r-t 

*41 


i 

MV 


valid  for  I 

Mmo* 


<U*) 
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Co^artng  both  equations  for©  end/.?1  ,  it  i.  sridsnt  that  th«  following 
condition  it  satisfied!  6 


^msobs«wic  ^ 

SUPfRSONlC 


)  -  @(MSubsow»c) 


<lUb) 

which  can  also  be  deduced  fro*  its  physical  significance.  In  addition  to 
this  behawier,  ^  has  a  non-rani  Aing  Halting  Talus,  if  the  Mach  nuaber 
approached  infinity,  which  will  be  of  importance  for  aany  future  calculations. 

■&tt*fe-)1=r  -  °-60264L,.«s) 


^  decrease  aonotonioally  from  1  at  Maoh  Ho.  1  to  the  limiting 

ralue  (I40)  at  Mach  Bo.  infinity.  The -  function  is  represented  in 
graphs  9#  10  or  11  by  the  curves  for  f  »  1. 

Comparing  )  11 )  and  (li*b)  one  obtains 

"  -H*u)  f°r  1  *  Mmav  (1U) 

Squat  ion  (l^ )  oombinsd  with  (U4)  and  (li|d)  results  ini 


fJiM . 

«  X 

7  ^(Mk) 


(15m) 


(15b) 


Bquation  (15b)  expresses  the  diffuser  effloienoy  of  the  subsonle  solution 
mathematically  by  the  'vf  function  at  the  supersonio  velocity  which  phyeioally 

*jd,t  lat*r)‘  Howerer,  this  squation  is  rsry  eonrsnisnt  for 
nuasrical  oaloulations* 
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Section  it  -  Discussion  of  the  Solutions* 


Complete  numerical  evaluation  of  the  solution  according  to  equations 
10,  12,  5,  5a,  15a,  15b,  2lu  is  represented  in  Fig,  2  to  S.  They  shew 
final  velocities  Ulet  ■,  Uip  (Fig.  2),  static  pressure  I»  and  total  pressure 
ratio  Kd  (Fig,  5,  k)$  total  pressure  ratio  Kjp  (Fig.  and  diffuser 
efficiency  'vjj  (Fig.  6,  7*  S)  of  the  ft  “solution  with  initial  velocity 
or  initial  Mach  number  Ma^  and  f  as  independent  variables  or  parameters. 

Fig,  2  represents  the  numerical  evaluation  of  equations  10  and  12, 

It  shows  toe  final  velocities  M204  and  M2pas  a  function  of  initial 
velocity  Ifc  with  parameter  f.  Also  the  total  pressure  ratio  Kp0 1  , 

K  is  given  as  a  function  of  and  f *  It  is  seen  that  for 
the  velocities  of  the  subsonic  solution  are  connected  with  the 

total  pressure  ratio  I  or  even  (see  also  Fig,  5)*  while  toe 

velocities  of  toe  supersonic  solution  H2O  1  are  connected  with  imaginary 
K  3>*  or  with  Kj**.  >  I  •  Because  this  latter  means  entropy  decrease 
(equation  1+b)  of  the  final  state,  the  supersonic  (  oi  )  solution  has  no 
physical  significance.  However,  before  discussing  this  problem  in  detail 
(Part  C),  the  following  paragraphs  (Parts  A  and  B)  deal  first  with  the 
discussion  from  the  mathematical  point  of  view. 


A,  Velocity  and  Total  Pressure  Ratio  of  the  0 t  -  and  fl  -Solution, 

1.  Some  Special  Cases, 

(a)  In  case  of  f*L,  which  is  equivalent  with  free  jet  length 
zero,  equation  (9)  is  reduced  tot 

For  in  the  supersonic  range,  this  gives* 

Map  =  ■gy  and  according  to  (12)  with  (li*b)  and  (lit)* 


M*«t  =  M, 


Kjot  ~  1  and  Kpp  -  ^jrj  >/mJ)  “  ^Ns  (M,J 


which  means,  we  have  either  no  change  of  state  or  a  normal  compression 
shock.  For  Hj'  in  the  subsonic  range,  it  results  int 


Mz4.  *  -jqy  j  MI(}  =  M,  and  correspondingly! 


1/  _  1 

"•*  «5t)'8(i)‘|WAl 


•nd  K,0  =■  I 


0* 
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This  means  we  have  either  no  change  of  state  or  a  normal  “rarefaction 
shock", (the  opposite  of  a  normal  compression  shock),  which  is  connected 
with  K  j>ci >  I  •  In  Fig.  2,  the  curves  for  f**l  form  the  exterior  boundaries 

of  a  lower  and  upper  realm,  which  contain  the  subsonic  and  supersonic 
solution. 


(b)  To  study  the  case  of  f-0,  which  corresponds  to  a  jet 
discharging  in  an  infinitely  large  reservoir,  we  first  multiply  equation 
(9a)  with  f  and  than  take  f-O,  which  yields  a  linear  equation 

("m]  y+i —  M')  =0 

Assuming  %<l^aa3C,  the  expression  within  parentheses  does  not  vanish  and 
we  obtain 

Mia  °0 

while  the  supersonic  solution  evidently  is  indeterminate  (see  Fig,  2), 

*  Hnflv  gives  no  determinate  value  for  M  .  For  determining  K2fl  a 
equation  (12)  fails,  but  equation  (6)  gives 

K-pp  -TfHiaiS.)  -  =o 

Because  P(K2(3  j  sx)  *  1,  it  follows 

Kpp  -P(H)*0  and  (17b) (17e) 

Discharge  of  a  free  jet  in  an  infinitely  large  reservoir  gives  the  total 
pressure  equal  to  the  static  pressure  of  the  jet,  a  result  well  known  for 
incompressible  fluids. 


(c)  In  case  of  Mai  mo°  ,  we  have  to  assume  f+O  and  insert 
in  equ  (10)  ®  M^ny  *  "/Hp  ,  which  gives  (see  Fig  2) 


*JWandM^ 


I 

Mm** 


(13a) 


The  total  pressure  ratio  is  given  by  equ  (12) s 

*  0  (13b), 

while  Kw  is  indeterminate.  The  diffuser  efficiency  Vpfb  is 
indeterminate  by  equ  (13),  but  can  be  calculated  from  equ  (15b),  which 
gives  the  limit* 


(18c) 

This  result,  as  important  as  it  is  simple,  shows  the  convenience  of  using 
our  definition  of  the  diffuser  efficiency.  The  advantage  in  practical 
application  for  higher  Mach  numbers  lies  in  the  fact  that  'pp  for  Mach 
numbers  above  3  and  for  f  values  >0.6  used  in  conventional  wind  tunnel 
design,  differs  only  5%  or  less  from  its  limiting  value  f  (see  Fig.  6). 
Evidently  (see  Fig.  3)  equ  (18c)  is  valid  for  f-0,  too.  Nevertheless, 
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the  point  lle<  • =»  ,  t -0  It  alngiiler,  ea  aean  in  the  dtscunlm  of  3?- 
(aeaPSaotlon\#  0,  i,  equ  (52)). 

2,  Relation  between  the  and  (3  «dlolutlon. 


We  investigate  mar  tha  relation  between  the  c<  -  and 
/3  -solution.  2qu  (11)  indicates  that  their  velocities  are  interconnected 
in  the  earns  way  as  the  velocities  before  and  behind  a  noraal  shock# 
Considering  the  total  pressure  ratios  K*ct.  and  1*0  according  to  (12),  and 
assisting  Miot  ^  %obx  md  using  (lUd)  and  (lit),  we  obtain* 


@0^1  =  a  is  f  ^  \ 


*  Kpp  or  plti*  )  po  p  (19*)#  (19b)# 

Th®  total  pressures  of  the  ®t  and  ft  '•solution  are  connected  by  the 
total  pressure  relation  of  a  noraal  shock  occurring  at  the  velocity 
M201  •  Iqu  (U)  with  (19b)  demon  strata  that  the  state  of  the  supersonic 
solution  is  connected  with  the  state  of  the  subsonic  solution  by  a  noraal 
0 oppress ion  shock  occurring  at  the  velocity  of  the  supersonic  solution# 

It  should  be  emphasised  that  this  statement  deals  with  the  relation  between 
two  possible  final  states#  The  relation  between  the  initial  and  either  one 
of  the  final  states  is  not  that  of  a  noraal  shock  (except  for  f^L),  as 
shown  in  the  following  paragraph# 

3.  Relation  between  the  Initial  and  Final  Velocity. 

In  order  to  give  a  simplified  representation,  we  exclude  from 
tills  paragraph  the  three  special  eases  dealt  with  in  paragraph  1,  namely 
j ML,  f*Q,  Mai  •  00  •  Then  we  may  prove  the  following  two  statements 
(a)  and  (b)i 

(a)  If  1  <  %  <  Mam  then 

i  <»•> 

In  other  words*  If  the  velocity  (1)  of  the  flree  jet  is  supersonic,  then 
the  velocity  after  the  transformation  process  of  the  supersonic  solution 
(2  ot )  is  greater  then  that  of  the  free  Jet,  while  the  velocity  of  the 
subsonic  solution  (2(3 )  is  loss  than  that  after  e  noraal  compression 
shock  in  the  jet. 


Rqn  (9)  can  be  written  for  the  solution  M**,  as  follows  1 

T>b +  M*  "  7  "Sr  +  “  7  M' 

The  right  side  is  larger  than  ,  as  oan  be  seen  by  transforming  into 

(j-<)ir  +  (^r-T^r-  ■)*><> 

* (f  -0>(f’ -0%  M‘Z  * wMflh  *■  w  *1 < 
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Hence  ^lcL  +Mz«  >75-,  +M,  and  because  -r+M  is  a  aonotonically 
increasing  function  for  M  >  1,  it  follows  M**)  and  according  to 
(11),  M2ft  <  -ji-  ,  where  is  the  velocity  after  a  normal  shock 
occurring  at 

(b)  If  0  <  Mx  <  1,  then 

~fhu  =  Mz(5<  M,<  '  <_Mr  <  MlaL  =  (20b). 

In  other  words*  If  the  velocity  (l)  of  the  free  jet  is  subsonic,  then 
the  velocity  after  the  transformation  process  of  the  subsonic  solution 
(2/3  )  is  smaller  than  that  of  the  free  jet,  while  the  velocity  of 
the  supersonic  solution  (2od)  is  greater  than  that  which  would  be  obtained 
after  a  normal  "rarefaction  shock"  in  the  jet  had  occurred.  The  first 
part  of  the  proof  for  the  solution  llzp,  is  the  same  as  for  M  za  above. 
^enco  TTift+  >  -rj  +Mi  9  but  because -r +M  is  a  monotonlcaUy 

decreasing  function  for  U  <1,  it  follows  &2{5  <  and  with  (li)i 
M  2ct  >  m:  .  Here  ^  is  the  supersonic  velocity,  from  which  the  velocity 
Hi  would  be  attained  by  a  normal  compression  shock.  In  a  mathematical 
nc*  physical  sense  we  can  say  also,  ~r,  is  the  velocity  after  a  normal 
"rarefaction  shock"  at  the  velocity  had  occurred. 

B»  Changes  of  State  in  the  Brea ^ure-Volo city  Diagram. 

1®  Representation  of  the  (3  -Solution® 


Fig.  3  is  a  dimensionless  pressure-velocity  diagram  and  gives 
a  representation  of  the  final  state  (M*p  ,  Pap  ),  which  belongs  to  each 
initial  state  (M^,  Pi)*  &ach  initial  state  lies  on  th©  upper  isen tropic 
line  s^  corresponding  to  a  total  pressure  ratio  K*L.  As  known,  the 
isen  tropic  lines  for  higher  entropy  eg  >  B'i  are  affine  to  the  upper  one 
by  being  multiplied  with  K<1,  according  to  equ  (5a).  We  designate  them 
by  the  K— value  which  is  identical  with  the  dimensionless  static  pressure 
P  at  zero  velocity  (M-0).  Only  the  subsonic  ( ft  )  solution  with  physical 
significance  Is  represented  here,  while  the  supersonic  ot  -solution  without 
physical  existence  is  omitted.  The  graph  contains  the  results  of  numerical 
evaluations  of  Jlap,  Kw  and  of  Pap  out  of  equations  (10),  (12),  and  (5a) 


All  of  the  final  states  of  the  (3  -solution  lie  in  a  region  which 
is  bounded  as  follows  s  By  the  P-exia  at  the  left-hand  sidsj  by  the  initial 
isentropic  line  from  M*0  to  M*1  at  the  upper  right j  by  the  well-known  "normal 
shock  polar"  from  K»1  to  M“  jcf'at  tlie  lower  right  j  by  the  M-axis  from  Um  tr 
to  M-0  at  the  bottom.  M™* 

The  final  states  (l Pg)  are  formed  by  the  intersection  of  two 
families  of  curves,  denoted  by  the  parameter  f  and  the  ip* -Hal  unr»y>  number 
Ma.  The  family  of  the  parameter  f  ranges  from  f-O.g  to  0.1*  The  curve  for 
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f«l  is  identical  with  the  above  mentioned  upper  right  border  (see  equ  (16b)) 
and  lower  right  border  of  the  region  (see  equ  (16a)),  while  the  curve  far 
f-0  is  identical  with  the  left  border  (17*),  excluding  the  origin  (fO, 

Ittu  »<*'  ).  All  f-curves  start  at  one  single  point  P’flL)  and  terminate 

(for  f+0)  in  one  single  point  (M«  ,  P-0)  according  to  equ  (13a),  (ISb). 

The  family  of  "initial  Mach  numbers"  ranges  in  thin  graph  from  0,2  to  U.5, 
in  order  to  obtain  clearness  in  representation,  while  evaluations  far 
initial  Mach  numbers  up  to  10  or  infinity  are  to  be  found  in  several  of  the 
following  graphs  of  (see  Fig.  6,  8)  and  other  wind  tunnel  characteristics. 
The  curves  for  "initial  Mach  numbers"  start  for  Ma^  4  1  at  the  upper 
isentropic  line  on  their  corresponding  isentropic  values  (l6b)j  far  >  1 
they  start  at  the  normal  shock  polar  on  their  corresponding  values  of  a 
normal  shock  (l6a).  The  case  when  f  approaches  one,  which  means  a  Vanishing 
afeert  jet  length,  is  therefore  identical  with  no  change  of  state  far 
Ma^  1  and  with  a  normal  shock  for  >  1.  The  M&j-curves  terminate 

at  the  Praxis  (equ  (l?a)  in  points  with  pressure  equal  to  the  initial 
static  pressure,  Pg  *  P^  (equ  17b)* 

2.  Graphical  Solution. 

The  pressure  velocity  diagram  offers  a  possibility  of  graphical 
determination  of  the  solution,  which  is  a  generalized  version  of  toe  well- 
known  method  for  a  normal  shock  process.  Due  to  the  fact  that  (~  a^js^onsf  =  w  § 

the  conservation  of  mass  and  momentum  can  be  represented  by  simple  graphical 
processes.  The  state  after  a  normal  shock  can  be  found  by  constructing  the 
tangent  to  the  initial  isentropic  line  at  the  point  of  initial  supersonic 
velocity  and  determining  the  second  point  of  contact  of  the  same  straight  line 
to  another  isentropic  line  in  the  subsonic  range. 


In  our  diagram  of  dimensionless  coordinates,  the  tangential 
relation  can  be  obtained  by  equ  (5)* 

H£iw. .  (a, 

Our  equation  system  (3)  with  f-a  describes  the  above  mentioned  graphical 
solution  for  a  normal  shock.  The  final  state  lies  on  the  normal  shock 
polar  with  X*X  NS  (%)• 

In  oase  of  f#L,  the  equations  (3)  state  that  toe  tangent  to 
the  isentropic  line  s^  at  the  point  of  initial  state  (Mj,  P^)  has  to  be 
rotated  counterclockwise,  until  its  slope  isjnultiplied  with  f  <  1,  and 
then  the  point  of  contact  of  this  rotated  straight  line  to  another 
isentropic  line  has  to  be  constructed.  This  new  point  of  contact  is  ths 
final  state  (Mg,  Pg),  The  designation  of  its  isentropic  line  gives  Kjj. 

It  can  be  seen  tram  the  graph  that  we  find  always  a  subsonic  solution  in 
the  above  described  region  of  Fig  3  ^ith  <  1  or  even  « 
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Additionally,  there  are  many  (M^,  f)  combinations  whore  a  second  point  of 
contact  at  the  right— 'land  side  of  in  the  supersonic  range  can  be 

found  which  lias  on  an  isen tropic  above  the  initial  isentropic  line*  This 
means  K  Po<  >  1*  In  many  other  (M^,  f)  combinations,  which  $.ve  M20(>Mraa3C 

(aee  Fig  2),  nojgraphical  representation  of  the  oc  -solution  is  possible, 
because  the  respective  pressure  is  imaglnsiy  (non-existent  gas)* 

3.  Approximation  for  Incompressible  Flow, 

Figures  2,  5j  6,  7,  3  show  that  tree  ft  ^solution  is  steady 
throughout  when  the  initial  jet  velocity  decreases  from  the  supersonic 
range  through  the  sonic  speed  down  to  subsonic  velocities  and  finally 
approaches  the  incompressible  region.  It  seer.®  desirable  to  compare  the 
evaluation  of  the  general  compressible  equation®  in  the  incompressible 
region  with  more  simple  formulas  which  describe  the  incompressible 
transforation  process  of  a  free  jet,  known  as  "sudden  widening  of  a 
pipe". 


The  incompressible  formulas  can  be  obtained  a a  follows* 
rhe  continuity acd  momentum  equ  (2)  become  for  ?2-  ?>  *? 

^5;  ” /•«'■?  >  -Sr=f  (22a)  (23b) 

Eirfnatlon  of  ^  gt™«,  p,  ,p,  . 

Introducing  the  total  pressure  of  the  initial  state 

„  i.  ?  w,1 

P°  ~  Pt  2  ~  results  in 

P»  *  t  ~2f( i-f))  (23) 

This  aquation  together  with  (22b)  presents  the  solution  (v*f  Pg)#  if  Wn, 
f  and  p0,  9  are  given.  If  we  prefer  dimensionless  values,  we  obtain, 
by  neglecting  terms  of  second  order,  the  approximate  solution* 

-f.-  I  — Ma,  W 

Denoting  the  total  pressure  of  the  final  state  with  pi  *  Pi. 

and  eliminating  wg  and  by  means  of  (22b)  and  (23),  we  obtain  the 

difference  in  total  pressures i 

P.-P.'-  (25) 

^or  tl»  total  pressure  ratio  we  obtain  trie  incompressible  approximation 

f)*T  Ma,1  ^  K>  -  I  -fl-f)1  MS  (26) 

^valuation  of  (26)  and  comparison  with  the  exact  values  of  (12)  give  for 
Mai  "  °*2  differences  in  Kp  between  3  x  10~5  at  f  *  0.3  and  5  x  lO-^  at  f  -  0 


Fig.  k  gives  the  enlarged  upper  left-hand  corner  of  the  pressure  velocity 
diagram,  shewing  i.he  exact  calculated  compressible  values  for  •  0.2 
and  0«l|  (same  as  in  ^ig.  3)*  compared  with  values  for  *  0.1,  0.2,  0«3| 
and  0j»,  obtained  by  the  inconqpressihle  approximate  formula  (2U)  above. 

The  comparison  is  an  excellent  check  for  the  results  of  the  general 
compressible  equations. 

C.  Proof  that  the  Diffuser  Efficiency  of  the  (i  -Solution  is  Smaller 
than  one  and  that  tl»  oi>  Solution  lias  no  physical  ^kiatence.  "** 

d_K  dH 

1,  Establishing  the  Equations  for  d±  and  A L _ i- 

Aa  already  explained  in  the  survey  of  this  Section  it,  it  can 
be  seen  by  numerical  evaluation  represented  in  Figs.  2,  5,  6,  7 ,  and  & 
that,  assuming  f  <  1,  K «,  is  either  larger  than  one  or  is  imaginary,  while 

always  Kp  <  1  or  even  %  <  1.  This  means  that  the  oC  -solution  has  no 
physical  significance.  However,  the  subsonic  ft  -solution  exists  physically 
and  shows  that  the  actual  diffuser  losses  are  larger  than  those  of  a  normal 
shock  at  test  section  velocity.  It  seems  desirable  to  prove  these  important 
results  directly  from  the  equations  in  addition  to  the  numerical  evaluations 
given. 


The  proof  will  be  made  using  the  fundamental  equations  (3), 
by  proving  first  that  >  l  and  <  |  •  Because  the  pressure  and 

mass  flow  functions  (equ  5)  are  transcendental  and  the  latter  cannot  be 
written  explicitly  in  the  velocity  M,  a  proof  with  these  special  functions 
is  very  inconvenient.  Instead  of  that,  a  mare  general  proof  will  be  given, 
using  a  generalised  form  of  isentropic  curves  in  a  pressure  (y)  versus 
velocity  (x)  diagram. 

In  the  following  only  positive  values  of  x  and  y  are 
considered.  The  equation  y^K  .  f  (x)  shall  represent  the  family  of  isen¬ 
tropic  lines  denoted  by  the  total  pressure  ratio  K,  with  F  (0)*<L.  F(x) 
shall  be  a  monot onically  decreasing  function  with  steady  first  and  second 
derivatives.  F'  (x)-0  is  allowed  for  x*0  only,  while  the  second  derivative 
shall  not  vanish  at  any  point.  Hence,  F(x)  can  be  concave  or  convex,  but 
any  inflection  point  is  excluded,  (x^,  y^)  is  an  arbitrarily  chosen  fixed 
point,  from  which  a  tangent  with  the  slope  a  is  drawn  to  one  isentropic 
line  out  of  its  family,  where  ytj)  is  the  point  of  contact.  We  define 
m  as  positive  for  clockwise  rotation  of  th«  tangent  around  (x^,  y^),  so 
that  a  >  0  covers  all  possible  cases.  (Xj,  y^)  does  not  necessarily  lie 

on  tee  initial  isentropic  with  K*l.  This  is  an  important  fact  for  later 
application  when  the  teat  chamber  pressure  is  not  equal  to  tee  Laval  nozssle 
exit  pressure.  Then  the  following  equations  are  valid* 


-y, 

xt~  xT 


=-m 


*,  yT*K-F(xT) 


j  — m  =  K*  Fj(Xt) 


(27a,  b,  c) 
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Equation  (27a)  corresponds  to  (3a),  (27b)  to  (5a),  and  (27c)  corresponds 
to  (Jb)  with  (5a).  We  consider  (x^,  y^)  as  fixed  and  m  as  independent 
variable,  while  (xy,  yy)  and  K  are  dependent  variables.  If  we  denote 
for  s  implif icat  i  on  that  =  9  (*)  and  g  as  the  inverse  function  of 

8*  (27o)  becomes  g(xT)  «  or  XT  }  and  (27b)  gives  yT  -  K- 

Elimination  of  xy  and  yy  in  (27a)  by  means  of  the  last  two  equations 
yields 

2  K-F[g(-gr)J -y,-mx,  +  =0 

Dividing  with  K  and  introducing  (u  =  -j-r  and  K  as  new  variables  with 


Ifli 

0  m 


gives  GfaK)  sFf§(( *) 

d 


3SL  , _ £ 

dK  K 


Now  we  obtain 


the  derivative 

dK  , 

(£).,_ 

1 

K 

dm 

(ffim 

-JZL  f 
KJ  ( 

1?)k  + 

(M) 

In  determining  (M.)  - -x,+(U  H. +5 

l  S(uyK  ^5  a<u  a 


we  have  to  consider  that  yg-  ~g(*r) 

abd  (U  -g(Xr)  ,  and  that  remains  finite  under  the  afore¬ 
mentioned  assumptions.  Hence,  it  be  comas  =>  §  —  /,  , 

Finally  w»  obtain 

fg  r. ,  -sy  y»  1  _  m  ft y>  i 

dK  K  f  “  K  i1+  m(*,-xr)  f  (20a) 

-^r  “  ~m  {'  “  }  ron<y  /»■  m*o)  (20b). 

^hese  aquations  shear  two  different  cases,  namely 

(a)  if  *1  >  Xy  or,  what  is  equivalent,  7i  <  7^,  then  -j~  >  0  (29a) 

(b)  if  <  Xy  or,  what  is  equivalent,  7^  >  then  <  0  (29b) 

If  we  now  insert  the  physical  quantities  of  our  case  and  consider 
that  ra  is  proportional  to  f  according  to  (3),  we  obtain  for  f*0 

~  ~f*{1  "%■}  **  (30*) (30b) 

or  with  the  special  value  of  m  according  to  (3)  in  (20a) 1 
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JL^±/I+  l 

d1  'M  f-®.  (*,-*>)/  (51). 

Determining  the  limit  of  (51)  for  f  approaching  zero,  we  have 
to  exclude  !»Sa^  *  00  and  may  then  use  (17a)  for  Mg  and  (17c)  for  P^,  which 


results  in* 


(Note  t  ^  Mmav  I ) 


ton  -4?- 


/*0 


(32) 


Taking  simultaneously  the  limits  of  f*0  and  Ma^  «<»  ,  the  quantity 
of  ^  ■  is  indeterminate*  It  might  be  mentioned  that,  violating  the 

assumptions  which  had  to  be  made,  (ISc)  would  give  for  f“0  the  value 
1,  while  (32)  would  give  for  Ma^  mco  the  value  1*0&j2S  (for  *y  ■  1*1|05)* 

Nevertheless,  this  behavior  is  without  practical  significance,  as  Fig,  fi 

shows* 


2*  Discussion  of  the  Equations, 

Investigating  the  sign  of  -§j-  or  4^f-ire  differentiate 

between  the  following  cases  for  the  oi-  and  for  the  3  solution  (f*0)i 

(a)  v If  is  in  the  supersonic  or  in  the  subsonic  region,  we  have  for 
f  <  1  always  >  U2p  according  to  (20a)  and  (20b),  which  results  according 
to  (29a)  in  >  0  and  ■  >  0,  Hence,  for  increasing  f,  Kt>e  and 

*7,9  must  increase,  too.  Because  for  f**l  we  have  r*(wal  ,  according  to 
equ  (l6e),  (16b),  it  follows 

for  f  <  1  becomes  ^  <  1,  (33) 

(b)  If  5^  is  in  the  supersonic  range  and  f  <1,  but  large  enough  to  give 
a  real  <*-  solution  with  Mict  <  (see  Fig.  2),  we  have  according  to 

(20a),  <  Uiai  ,  which  results  according  to  (29b)  in  -^Tjp0  K  0# 

Hence,  for  increasing  f  from  values  f  <1  to  f"l,  Ka*  must  decrease 
until  "1  is  reached  far  f"GL  (equ  (16a))  and  it  follows 

far  f  <  1  becomes  Kj*  >  1  (34») 

The  case  of  MaoC>  Mp^-r  involves  imaginary  pressures  and  is  of  no  physical 
interest* 

(c)  If  IL  is  in  the  subsonic  range,  we  also  need  to  consider  only  solutions 
with  Mi(1(.<fima3C,  In  this  case,  it  becomes  lf^<M20(  far  f  <  1  (according  to 
20b))  and  also  far  f»l  (according  to  16b)),  which  results  in-gj-<0 
(according  to  29b)*  Hence,  with  f  decreasing  from  f*4L,  it  follows  that 
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K  ot  increases  from  its  value -^->1  for  fl  (see  16b)),  which  moans 
far  f  4  1  becomes  Kp*  >  1,  (34b) 

Evidently,  in  both  cases  (b)  and  (c),  if  is  fixed,  Kj,  increases 

raonotonically  with  decreasing  f  until  Kpod  »°°  is  reached  (see  Fig.  2), 
Decreasing  f  further  gives  M2t(>  which  results  in  imaginary  pressure 
and  imaginary  total  pressure  ratio. 

(d)  0^  interest  may  be  the  special  case  of  Ma^  ma°  0  which  involves  the 

(se®  eclu  (S))  and  therefore  transforms  equ  (31)  for  f*0  into 
j'l  “  v  •  integration  with  the  arbitrary  constant  chosen  so  that  'V-ofi  "1 
for  f"l,  results  in  "(vp  *£  in  accordance  to  (lgc).  1 

Equations  (33)  and  (34a),  (34b)  prove  the  statements  made  in  the 
beginning  of  this  Section  3  G. 

Now,  we  can  differentiate  between  the  three  regions  of  solutions 
in  the  Mg  va«  diagram  (Fig.  2),  exduding  f»l« 

(a)  0  4  M2Ij<1  are  solutions  with  physical  existence,  because  their 
change  of  state  is  connected  with  entropy  Increase,  due  to  K-nn  <  1 
or  even  1*5  fld, 

(b)  1  <  are  solutions  without  physical  existence,  because  they 

belong  to  a  change  of  state  with  Kj#  >1,  or  entropy  decrease,  violating 
the  second  theorem  of  thermodynamics.  The  final  state  Itself,  however, 
belongs  to  a  gas  with  physical  existence. 

^  "max  ^  M2ot  are  solutions  without  physical  existence,  because  the 
final  state  belongs  to  a  non-existent  gas  with  imaginary  pressure. 

Therefor®,  only  the  subsonic  solution  has  physical  existence. 

Or,  in  other  words*  If  the  initial  velocity  in  the  Jet  is  either  subsonic 
tsr  supersonic,  the  final  velocity  in  the  diffuser  after  completion  of  the 
transformation  process  is  subsonic  in  all  cases. 

This  careful  study  fltthcsrning  the  physical  existence  was  necessary 
because  it  cannot  be  stated  in  advance  that  the  supersonic  solution  will 
have  no  physical  existence.  For  instance,  we  will  se©  in  another  report 
(Ref.  26)  that  in  certain  cases  where  the  test  chamber  pressure  is  not  equal 
to  nozalo  exit  pressure,  the  supersonic  solution  exists  physically  in  addi¬ 
tion  to  the  subsonic  one.  ^ 

3.  Conclusions*  The  Physical  Nature  of  the  Diffuser  Process. 

Concerning  the  diffuser  flow  process  the  following  conclusions 
can  now  be  drawn.  The  final  velocity  after  the  transformation  is 
than  that  attained  after  a  normal  shock  (equ  20a).  Equ  (33)  states  that 
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the  diffuser  losses  are  always  larger  than  those  of  a  normal  shock  at 
Jet  velocity*  Jfence,  the  diffuser  process  is  not  that  of  a  normal  shock 
at  Jet  velocity.  Neither  does  it  correspond  to  that  of  a  normal  shock  at 
a  higher  velocity*  to  which  the  jet  velocity  is  accelerated  by  iaentropic 
expansion  from  F.,  to  a  larger  area.  This  was  once  assumed  by  II,  Ram 
(Ref,  2  and  5)  for  the  purpose  of  discussing  his  experimental  data  on  the 
pressure  distribution  along  the  diffuser  and  proved  by  him  to  be  false* 

In  the  pressure  velocity  diagram  (*’ig.  3)  the  final  states  lie  definitely 
on  one  side  of  and  not  on  tire  normal  shock  polar  (Sec,  1*,  B,  1), 

Neither  is  the  process  an  addition  of  losses  in  tto  jet  or  on 
its  surface  followed  by  a  normal  shock*  as  often  presumed  for  explanation 
of  the  low  diffuser  efficiency.  If  the  second  process  is  a  normal  shock, 
the  state  upstream  of  it  can  be  only  the  <*-  solution  (see  Sec,  I4 ,  A,  2), 
which  has  K*  >  1,  or  total  pressure  increase.  Hence,  instead  of  jet 
processes  with  total  pressure  losses,  such  ones  with  total  pressure  gain 
(entropy  decrease)  should  precede  a  normal  shock,which  leads  to  the  final 
state  calculated  here.  This  same  conclusion  is  valid  also  for  the 
assumption  of  any  combination  of  oblique  shocks  followed  by  a  normal  shook 

Neither  is  the  process  an  addition  of  a  normal  stock  at  jet 
velocity,  followed  by  a  recapture  and  transformation  process  of  a  free 
jet  with  subsonic  velocities.  Such  process  would  lead  to  higher  total 
pressure  ratios,  because  to  smaller 

velocities  and  higher  static  pressures,  as  can  be  seen  in  Fig,  3, 

It  seems  useless  to  search  for  explanations  of  this  diffuser 
flow  process  after  a  jet  expansion  by  means  of  a  normal  shock.  The 
process  here  investigated  is  a  more  general  process  between  the  areas 
and  )  ^1*  from  which  the  normal  shock  is  only  a  special  case  for 

F2  approaching  F^.  The  special  ease  never  can  explain  the  general  case. 
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Section  5*  Calculation  of  Other  Wind  Tunnel  Characteristics , 

A.  Dimensions  of  the  Diffuser  Thirst  and  Teat  Chamber  Pressure  Control, 


The  dimension  of  the  second  throat  has  been  always  an  important 
feature  in  the  design  of  supersonic  wind  tunnels.  The  result  obtained 
above  (Sec.  1*,  C,  2),  that  the  velocity  in  the  diffuser  entrance  section 
is  subsonic  after  completion  of  the  transformation  process  is  very 
important  to  the  following  discussion.  Assuming  a  certain  restriction 
of  area  in  the  central  part  of  the  diffuser  and  a  pressure  maintained 
below  a  certain  limit  downstream  of  the  diffuser  by  compressors  or  other 
means  (as  a  vacuum  tank),  it  fellows  that  the  diffuser  acts  as  a  conven¬ 
tional  Laval  nozzle  with  accelerated  flow  throughout .  Therefore,  the 
diffuser  has  subsonic  velocity  in  its  convergent  part,  sonic  velocity 
in  the  minimum  area  and  supersonic  velocity  in  its  divergent  part. 

Because  the  minimum  area  of  the  diffuser  F*‘  acts  as  a  second  throat 
following  the  first  throat  F*  of  the  Laval  nozzle,  its  size  is  determined 
by  F#  and  the  total  pressure  losses  in  between,  the  equations  of  continuity 
and  energy  give* 

P*!  P.'  « 

F*>“  9*  “  p*  Po  (35) 


The  total  pressure  ratio  involved  here  is  according  to  (12),  because 
the  flow  in  the  convergent  section  of  the  diffuser  downstream  of  the 
transformation  zone  is  an  accelerated  subsonic  flow  which,  as  known,  has 
practically  no  losses. 


Using  ®,  (see  (2b))  and  ^  according  to  (13),  we 


transform  (35)  into 

Ma,>l  and  (36) 

©  cr  $  are  functions  of  M,  according  to  (5)  or  (ll+)j  ^  is  a 
function  of  f  and  L5^  according  to  (13)  or  (15a),  (15b).  Sometimes  it  is 
convenient  to  use  another  ratio,  deduced  from  (56)  j 


F*> 


for  Ma>l 


and 


h  *  > 


for  Mat  <*  I 


(37) 


Fig.  9  shows  equ  (36)  as  a  function  of  f  and  nozzle  Mach  No.  from  0  to  5, 
inclusive  of  the  limit  for  Mach  Ho.  infinity*  Fig,  10  shows  equ  (36)  as 
a  function  of  f  and  noftfl*  velocity  in  the  reduced  scale,  to  enable 
representation  of  small  f.  Fig,  13.  shows  equ  (37)  as  a  function  of  f  and 
nozzle  Mach  No.  from  0  to  5*  ^h©  limit  for  Mach  No,  infinity  is  the  same 
for  all  f-valuee  and  is  equal  to  the  limit  of  (equ  (lUc).  The  curves 
for  f*l  In  Fig,  9#  10,  11  represent  the  -function,  due  to  ^  *1. 
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If  the  f-value  and  Mach  No.  in  the  test  section  are  given, 
only  one  specific  dimension  of  the  second  throat  is  determined  by 
the  equation  (36)  to  which  quantity  the  adjustable  diffuser  throat 
has  to  conform,  in  order  to  fulfill  the  fundamental  assumption 
made  in  this  analysis,  namely:  Equilibrium  between  the  pressure  of 
test  chamber  and  Laval  nozzle  exit.  Analysis  in  another  report 
(Ref.  26)  of  the  case  p  ^P]_  shows  that  in  a  certain  range  Pc<Pi 
results  also  in  a  subsonic  solution  but  with  Kr  smaller  than  for 
pressure  equilibrium,  while  PC7P1  results  in  %  larger  than 
for  equilibrium.  This  fact  causes,  according  to  (35)#  an  F**  larger 
and  smaller,  respectively,  than  that  for  pressure  equilibrium.  Or, 
expressed  the  other  way:  Varying  the  adjustable  diffuser  "throat, 
changes  the  test  chamber  pressure*  One  specific  dimension  of  the 
throat  according  to  (36)  gives  chanber  pressure  equilibrium  with  the 
nozzle  exit  pressure,  while  an  enlarged  throat  gives  smaller  chamber 
pressure,  and  a  reduced  throat  gives  higher  chamber  pressure  than 
the  nozzle  exit  pressure.  It  may  be  emphasized  here  that  the  control 
of  the  test  chanber  pressure  by  means  of  the  adjustable  diffuser  as 
outlined  above  is  an  exact  result  of  a  one-dimensional  calculation 
without  consideration  of  friction.  No  assumption  of  pressure  transfer 
upstream  through  subsonic  boundary  layers  is  necessary,  as  is  sometimes 
done. 

Evidently  all  conditions  discussed  above  are  independent  of  the 
downstream  pressure  as  long  as  it  is  below  a  certain  limit,  so  that 
sonic  velocity  in  the  minimum  cross  section  with  maximum  mass  flow 
per  unit  area  is  just  maintained.  Therefore,  pressure  equilibrium 
in  the  test  chamber  equivalent  to  a  parallel  jet,  which  is  very  im¬ 
portant  for  actual  wind  tunnel  testing,  is  independent  of  fluctuating 
intake  pressure  at  the  blowers  for  continuous  operating  tunnels.  It 
is  also  independent  of  vacuum  tank  pressure  increasing  with  time  for 
intermittent  suction  type  tunnels,  ^o  guarantee  safe  operation, 
actually  an  acceleration  to  a  slightly  supersonic  flow  is  required  at 
least  for  a  short  distance  downstream  of  the  second  throat,  with  trans¬ 
ition  to  final  subsonic  velocity  hy  a  normal  compression  shock  at  a 
low  supersonic  Mach  No.  with  very  little  total  pressure  loss. 

B.  Volume  Flow  into  Blowers. 

The  volume  flow  V  per  unit  time  into  the  blowers  is  an 
important  characteristic  in  the  power  plant  layout  of  supersonic 
wind  tunnels.  It  can  be  calculated  directly  from  the  diffuser  effi¬ 
ciency  under  the  assumption  that  all  additional  losses  occurring  between 
second  throat  of  the  diffuser  and  the  blower  intake  are  negligible 
compared  with  the  losses  which  have  already  occurred  at  the  beginning 
of  the  diffuser.  This  is  fairly  valid  far  a  reasonable  subsonic 
diffuser  and  straight  piping  without  corners  or  coolers  far  supersonic 
test  section  velocity,  as  measured  in  the  Peenemuende  wind  tunnel 
(Ref,  3).  Further,  the  intake  velocity  into  the  blowers  is  considered 
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to  be  ao  small  that  the  static  pressure  and  temperature  can  be  replaced 
by  their  stagnation  values  p#  and  T#,  with-^*  - 

Xn  the  following  we  make  the  special  assumption  that  no 
intercooling  shall  be  provided  between  test  section  and  blowers* 

Hence  the  compressor  inlet  temperature  is  equal  to  the  stagnation 
temperature  of  the  air  at  the  exit  of  the  fixed  diffuser  a  nd  to  that 
before  entering  the  test  section,  or  T#  •  T  .  In  case  intercooling 
is  provided,  the  formulas  of  paragraphs  B,  C,  and  D  liave  to  be 
changed  accordingly* 


which  gives  with 


The  continuity  of  mass  flow  is 
JW,  =  F*f  a* 

Si  Pe 


expressed  by 


and  with  equation  (lip) 


JL  _  9*  „*  & 

5  ?o  “fc 

C8? 

J*  _  (  2  y  «->  fr  „ 

~  v?+r7'  ifa  a° 
velocity  of  sound  &0  is  taken  at  the  stagnation  temperature  TQ» 

^  be  taken  from  Fig.  9  according  to  equ  (36).  The  function  4> 
(36a)  has  to  be  replaced  by  ©  in  case  of  subsonic  velocities* 


For  rxusssricaX  evaluation  we  choose  T  *  295*  K*  and 
&e  in  (n/ses)  by  \  ■  331.3  VJf  ,  which  corresponds  to \  in 

and  in  (h?).  With  Y  *  1*U05  the  result  is 


V«  f 

'  "Awl 

|  =*  198.8- 

F,  i 

J 

*1* 

(36b) 

Figs*  12  and  13  are  plots  of  (36b),  a a  a  function  of  Mach  No.  and  with 
parameter  fj  Fig.  12  for  Mach  No.  up  to  5,  Fig.  13  in  a  ami-logarithmic 
scale  up  to  8*5  'with  the  llxiit  tot  Mach  No*  infinity*  Per  revalues  of 
conventional  wind  tunnel  design  (f  between  OS  and  0.5),  the  V#«eurves 
above  M&*1  first  show  a  slight  increase  with  a  maximum  between  Mach 
No.  1  and  2,  whose  location  depends  on  the  f- value  ,  then  a  slight 
decrease  and  soon  the  curves  approach  the  limit,  The  fact  that  V# 
changes  only  slightly  in  the  whole  supersonic  range  is  of  great 
importance  for  engineering  purposes  in  the  layout  of  the  blowers. 
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For  instance,  for  f-O.&i  we  obtain  the  following  -waliee* 
at  Ma-1*  n=f  -  224— ~ic  j  at  Ua*L*li225i  -§•  =  235  (suudaaa)j 

at  Ma*«>i  =  >87  ~£(c 


C,  Energy  Consumption  of  the  Bltwcers 

The  energy  consumption  of  the  blowers  can  be  calculated  by 
the  compression  ratio  of  the  blowers  -Es  =  -1-,  either  for  adiabatic  or 

isothermic  compression.  *s  long  as  one  blower  unit  may  be  with  several 
internal  stages,  can  be  used  ( approximately  up  to  Mach  No*  2.5),  the 
adiabatic  compression  is  applicable*  If  for  higher  Mach  numbers  two  or 
more  blower  units  in  series  have  to  be  applied,  intercooling  is  necessary* 
Hereby  the  compression  process  approaches  mere  nearly  an  isothermic 
compression  as  a  greater  number  of  external  compresses*  stages  are  installed, 


We  start  from  the  well-known  equation  of  compression  work  per 
unit  mass  for  a  complete  cycle,  with  the  same  assumptions  as  in  Section 
B  (compressor  inlet  temperature  T#  “  TQ) 

WiS=  R  T0-  In  and  Wad  =  -JL-  ‘R-T<s 


[(£)  - > ] 


We  multiply  by  the  factor  of  mass  flow  per  unit  time  and  unit  nozsle 
exit  area,  -Vy--^  and  obtain  the  compression  energy  per  unit  noaale  exit 


area. 


Lis 

F, 


(»») 


and  correspondingly. 

We  choose  the  sound  velocity  aQ  for  the  stagnation  temperature 
T0  •  293°K  and  the  stagnation  pressure  P0  -  1  atm  »  10332  and 


we  convert  into  international  kilowatts  by  the  relation  1  int  W  • 

102*00  -°y^-  .  We  finally  obtain  the  blower  energy  in  int.  KW  per  s£, 
assuming  100%  blower  efficiency* 


Lis  fint.Kw] 

F.  L  m*  J 


20140  ©■  In 


Kp 


and  -^d=  20140 — 


r  a_  -i 


(39b) 


Fig.  31*  shows  equ  (39b)  as  a  function  of  Mach  number  with 
f  as  the  parameter*  For  incompressible  flew  (Mach  number  approaching 
aero),  isothermic  and  adiabatic  energy  are  identical*  Fig*  15  showB 
the  curves  for  f"0*&*  in  enlarged  scale.  For  conventional  wind  tunnel 
design  (f  between  OB  and  0*5)  the  energy  curves  have  a  peak  between  Uadi 
No,  2  and  3,  Fig,  li*  shows  also  for  the  purpose  of  comparison  the  curve 
for  the  kinetic  energy  l^in  in  the  test  section  for  the  same  conditions* 

It  might  be  noticed  that  for  larger  f  at  subsonic  and  small  supersonic , 
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Mach  numbers  is  larger  than  or  L^,  but  for  higher  supersonic 

Mach  numbers  the  opposite  is  valid#  The  blower  energy  at  no  pressure 
recovery  (f«0)  is  equal  to  the  kinetic  energy  of  the  jet  in  the 
incompressible  range  only  (see  equ  (l*2d)),  while  in  the  compressible 
range  the  blower  energy  is  larger# 


The  kinetic  energy  is  calculated  by  U.v,  F*$*a*w2 


which  gives 


L  Kin 

F, 


(M-tr) 


r+» 


a0p0®M4 


(i*Oa) 


and  with  the  same  conditions  for  pQ,  TQ  as  chosen  above* 

“TT  [^1  -  TT  2OMO'0 -M* -  II770  ©M!  (tob) 


The  last  numerical  factor  is  the  kinetic  energy  for  Mach  No.  1#  Per 
Ma*  VT  the  maximum  value  «  19320  | is  attained# 


D#  Energy  Ratio 

As  explained  before  (Sec#  1,  C,  J)  the  energy  ratio  depends 
on  a  specific  but  arbitrarily  chosen  compression  and  cooling  process# 

We  use  the  same  assumptions  regarding  the  cycle  as  outlined  in  Section 
5,  B.  Hcwevera  the  values  of  these  ratios  are  independent  of  the  special 
choice  of  pe,  ”Q.  According  to  an  isothermic  or  adiabatic  compression 

process*  we  have  to  differentiate  between  the  isothermic  and  adiabatic 
energy  ratio.  Considerable  differences  exist  between  those  two  terras  in 
trie  higher  supersonic  range# 


Dividing  (!*0a)  by  (39a)  we  obtain  the  ieethermic  energy  ratio 

C  - 

^Ns 


L;< 


tr 

er+i 


M 


wd  the  adiabatic  energy  ratio 

c 


•ad 


Lkin 


M* 


ftF-J 


(la*) 


tfllb) 


16  is  a  send-logarithadc  plot  of  equ  (Ida)  and  (l^b)  as  a 
function  of  Mach  number  with  parameter  f  up  to  Mach  No.  10.  There  is  no 
finite  limit  for  Mach  No,  infinity.  The  energy  ratio  is  infinite  for  f*l 
in  the  whole  subsonic  rang®,  because  this  corresponds  to  total  pressure 
ratio  1  (equ  (lob) )  with  ho  compression  work  needed.  In  the  supersonic 
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rang®,  however.  represents  the  normal  shock  (equ  (16a))  with  finite 
losses  and  finite  energy  ratio.  When  f-O,  which  means  no  pressure 
recovery,  we  obtain  £  •!  for  Mach  number  approaching  aero,  but  <  1 
at  compressible  flow.  It  can  be  seen  that  In  the  supersonic  range  tbs 
C,  -values  decrease  rapidly  and  finally  all  become  less  than  one. 

Mr  conventional  wind  tunnel  design  corresponding  to  t  between  0.5  and 
0.8,  this  happens  between  Mach  number  5  and  5  for  both  the  adiabatic 
and  isothermic  energy  ratio. 

It  might  be  mentioned  that  Ackeret’s  definition  of  the 
energy  efficiency  of  the  diffuser  (Sec.  1,  C,  1)  can  be  written  as 


,  hence  becoming  negative  in  a  broad  supersonic 


range  where  1,  oven  though  pressure  recovery  is  actually  present 
(i  >o)  .  therefore  the  value  of  this  definition,  given  primarily  for 
the  compressible  subsonic  region,  seems  doubtful  in  tbs  supersonic 


range. 


In  the  subsonic  range  we  must  note  that  the  transformation 
losses  in  the  entrance  of  the  diffuser,  which  are  represented  in  fig. 

16,  are  beooming  comparatively  small  and  the  losses  in  the  subsonic 
diffuser  itself  and  in  the  other  parts  of  the  circuit  (comers,  honeycomb) 
which  are  not  included  here,  cannot  be  neglected  as  in  the  supersonic 
range.  Therefore,  experimental  values  of  energy  ratio,  which  include 
all  such  losses,  must  be  considerably  smaller  than  the  calculated  values 
in  fig.  16.  The  highest  known  energy  ratio  of  an  incompressible  free 
jet  wind  tunnel  is  that  of  the  California  Institute  of  Technology  in 
its  former  arrangement  with  a  very  short  free  jet  length  of  jr  "O.8. 

Clark  invHVan  (Ref.  12)  reported  an  overall  energy  ratio  of  5*6, 
referred  to  motor  output,  which  will  correspond  to  approximately 
£  <*k.8,  referred  to  blower  output.  The  calculated  energy  ratio 
based  on  the  transformation  losses  of  this  short  jet  (f*  approx. 

0.80,  Ma-0.23)  la  C  m  2k,  according  to  fig.  16 

In  the  case  of  incompressible  flow,  Ma-*0,  we  have 

Li's  =  L«d  ®  w(po-Pe)  (i*2a) 


and  —  g  §o  w 


which  gives 


j9oW* 

Pc-P* 


(i*2b) 

(1*20) 


Equation  (l*2a)  can  be  obtained  also  from  (59b)  as  well  as  equation  (l*2c) 
can  be  obtained  frora(Ula)  by  a  limiting  process  Ma  -*■  0,  considering 

that  In  -It  -*•  — °-~  p?  .  Using  the  difference  in  total  pressures 

»  ^  Po 
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expressed  by  f  (equ  (25)),  the  energy  ratio  (i|2c)  for  incompressible 
fleer  becomes  finally 


S  ~  Cad  * 


(i*2d) 


Now,  it  will  be  analyzed  how  the  needed  blower  energy  depends  upon  the 
length  of  the  free  Jet  test  section.  For  this  purpose.  Fig.  17  shows 
the  reciprocal  isothermic  energy  ratio  -tif  according  to  (IfLa)  and 

(U2d)  as  a  function  of  the  relative  free  jet  length  with  the  Mach  No, 
as  parameter  varying  between  0  and  10.  In  this  graph,  the  parameter 
f  used  in  all  previous  plots  is  eliminated  and  4r-  is  substituted 

with  the  aid  of  equ  (1*3),  using  for  the  expansion  angle  £  an 
experimental  value  of  e  -0.06l  (see  Sec.  6,  A). 


The  subsonic  curves  start  at  zero  energy  for  free  jet 
lengthy  zero.  For  incompressible  flew  (Ma«0)  the  blower  energy  for 
mail  jf  increases  approximately  parabolic,  according  to  (l+2d)  and 
(U3)«  The  supersonic  curves  start  at  a  finite  blower  energy  for  free 
jet  length  zero  to  overcome  the  normal  shock  losses,  and  increase 
nearly  linear  with  the  jet  length.  Hence  it  follows  that  for  sub” 
sonic  and  especially  low  speed  wind  tunnels,  increasing  the  free 
jet  length  increases  the  required  blower  energy  considerably.  For 
supersonic  velocities,  however,  and  especially  for  high  Mach  numbers, 
increasing  the  free  jet  length  increases  the  blower  energy  by  only  a 
relatively  small  amount.  For  instance,  increasing  the  free  jet 
length  from  1  to  2  increases  the  blower  energy  at  Ma«0  approximately 
200%,  at  Ma»l  approximately  50%,  at  Ma4u5  only  3%,  and  at  Ma-10 
only  5%.  Reducing  the  length  of  free  jet  in  order  to  keep  the  blower 
energy  small  is  therefore  important  for  subsonic  wind  tunnels  only. 

In  supersonic  wind  tunnels,  especially  for  hi^i  Mach  nunbers,  the 
length  of  the  free  jet  need  not  be  considered  from  the  standpoint 
of  saving  power  but  can  be  chosen  in  favor  of  optimum  testing  con¬ 
veniences  . 
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Station  $.  Discussion  of  Experimental  Results 

A.  Experimental  Values  of  Art*  Ratio  f 

A a  seen  in  the  preceding  stations,  the  tree  ratio  f  is  the 
important  parameter  of  all  the  theoretical  calculations  and  disoussions. 
Aotually,  it  represents  the  boundaries  of  our  one-dimensional  problem, 
the  nosele  exit  and  the  diffuser  intake  oross  seotioi^”  If  f  is  given  in 
a  speoial  case,  the  flow  prooess  through  the  diffuser  and  its  effieienoy 
dib  be  oaloulated.  Another  problem  is  which  f— 1 value  belongs  to  a  certain 
relative  length  1_  of  the  free  jet.  Or  expressed  more  exactlyi  Ifhat 

H 

diffuser  intake  area  is  necessary  to  produoe  a  smooth  recapture  process 
of  the  jet  of  a  given  length  at  pressure  equilibrium  in  the  test  ohamber? 

The  answer  to  this  question  shall  be  given  in  the  present  report  purely 
experimentally  and  is  represented  in  Fig.  18  as  a  plot  of  f  vs  1_  from  data 

used  in  different  German  free  jet  wind  tunnels,  which  have  been  in 
successful  operation  and  testing.  In  order  to  replaoe  the  relation  between 
f  and  1  by  one  value,  the  expansion  angle  £  is  introduced  by  the  equation 
ff  (see  Fig.  1) 

^  ~  (\  +  2  C-&-)2  (^3) 

This  formula  assumes  similarity  of  the  jet  expansion  in  reference  to  the 
trailing  edge  of  the  Laval  nozzle,  as  known  by  the  two-dimensional  analysis 
of  the  jet  boundary  by  W.  Tollmien  (Ref.  8), 

The  following  explanations  are  given  for  the  data  shown  in  Fig.  IS. 

The  Goettingen  AVA  supersonio  free  jet  wind  tunnel  11  cm  x  13  cm  was 
operated  with  model  and  bell  mouth  at  an  expansion  angle  £ •  0.053  ty 
Walchner  (Ref.  13).  The  Goettingen  KWI  supersonio  wind  tunnel  6  om  x  S  om 
was  equipped  with  recapturing  knife  edges.  The  test  section  was  half  closed  - 
half  open.  Without  model  the  expansion  angles  for  the  free  jet  boundaries 
was  determined  by  K.  Oswatitsch  (Ref.  li*)  to£*  0.06l.  For  the  purpose  of 
comparison  in  Fig.  IS  f  is  calculated  with  this  e-value  from  equation  (i+3) 
for  a  completely  free  jet. 

The  Aaohen  AIA  free  jet  supersonio  wind  tunnel  10  cm  x  10  cm  (Ref.  15) 
was  operated  by  the  author  with  model  and  bell  mouth  at  an  expansion  angle 
of  £  -  0.06l.  The  Peenemuende,  later  Kochel,  free  jet  supersonic  wind  tunnels 
I  and  II  of  I4O  om  x  UO  cm  were  designed  by  the  author  (Ref.  16)  for  long 
models  with  a  free  jet  length  of  1/11  •  1.75  end  with  a  bell  mouth  of  6  ■  0.071. 
The  Maoh  number  range  is  1  to  3*2  for  the  Goettingen  and  Aaohen  tunnels  and 
up  to  U*U  for  the  Peenemuende  and  Koohel  tunnels.  The  fact  that  all  data 
in  Fig.  IS  lies  close  together  proves  that  £  is  fairly  independent  of 
different  variables  such  as  bell  mouth,  knife  edge,  model,  Reynolds  number 
and  Maoh  number.  Even  the  experimental  range  of  several  low  speed  subsonic 
Goettingen  tunnels  with  free  jet  with  £  varying  from  0.033  to  0.066  aooording 
to  L.  Prandtl  (Ref.  17)  is  close  to  the  supersonio  data. 


AFTR-635U 


35 


In  addition  to  the  test  points,  two  curves  are  drawn  in  Fig.  lg 
according  to  equation  (43)  with  the  values  of  t  »  O.O&L  and  £  «  0.071 
which  cover  most  of  the  supersonic  data.  A  theoretical  understanding 
of  this  experimental  value  of  £  can  be  obtained  with  the  aid  of  velocity 
distribution  in  the  turbulent  mixing  zones.  The  curves  can  be  used  to 
determine  f —  values  for  free  jet  length  l/H  values  other  than  those  ex¬ 
perimentally  established.  Because  (I4.3 )  is  based  on  two  dimensional  flow, 
these  curves  in  Fig.  IS  will  become  gradually  invalid  for  larger  l/H  values. 
They  might  be  used  up  to  approximately  l/H  ■  4,  according  to  the  velocity 
distribution  calculated  hy  A.  Keuthe  (Ref.  9;  its  Fig.  9,  n).  This  also 
is  approximately  the  region  where  the  interior  core  of  undisturbed  velocity 
disappears  according  to  experiments  of  P.  Ruden  (Ref.  IS)  and  A,  Kuethe 
(Ref.  9). 

For  conventional  wind  tunnel  purpose  a  relative  free  jet  length  of  2 
or  3  seems  to  be  the  upper  limit.  For  the  purpose  of  application  of  this 
theory  to  turbo-jet  or  ram-jet  test  stands,  where  sometimes  longer  jets  are 
required,  the  curves  might  be  extended  with  caution  into  the  dotted  regions. 

It  shall  be  mentioned  that  the  increase  of  £  from  0»06l  in  the  Aachen  design 
0-934)  to  0.071  in  the  Peenemuende  design  (1937)  was  made  by  the  author 
intentionally  to  overcome  some  difficulties  in  establishing  pressure  equi¬ 
librium  at  Mach  number  around  1.5  when  large  models  were  installed.  In 
all  other  cases,  the  curve  for  £■  0.061  can  be  used. 

B.  Pressure  Distribution  Measurements  along  the  Diffuser 

As  explained  before  (Sec.  1,  B)  the  present  analysis  was  instigated 
by  the  results  of  the  pressure  distribution  measurements  along  the  diffuser 
of  the  Peenemuende  wind  tunnel  by  H.  Ramn  (Ref.  2)  which  now  shall  be  com¬ 
pared  with  the  theoretical  results.  The  author  himself  designates  in  Part  I 
these  measurement  results  as  "preliminary" 5  "the  purpose  of  the  investigation 
***  was  chiefly  the  testing  of  small  pressure  gages j"  "Since  it  was  realized 
that  measurement  techniques  were  still  far  from  perfect  ***  the  results 
should  be  used  only  as  a  general  survey." 

Therefore,  a  quantitative  agreement  between  theoretical  and  experimental 
results  cannot  be  expected.  However,  quoting  his  major  qualitative  results 
will  show  an  interesting  and  complete  confirmation  of  our  theoretical  con¬ 
clusions.  We  shall  disregard  all  experimental  results  dealing  with  the  so- 
called  "throttling  flap",  and  omit  those  concerning  various  test  chamber 
pressure®.  In  Ramm’s  recond  report  (Ref.  3),  Sec.  No.  4,  Experimental  Results, 
he  says,  "The  figures  (3,  4}  5  &  6)  show  the  following  characteristics  valid 
for  free  jet  wind  tunnel:  A  pure  subsonic  flow  in  the  diffuser  never  occurs 
at  the  test  Mach  numbers  1.22,  1.56,  1.S6,  and  2.50.  Transition  to  a  pure 
supersonic  flow  takes  place  only  with  a  diffuser  opening  which  is  very  large 
compared  with  the  free  jet  cross  section.  In  all  other  cases,  the  normal 
flow  pattern  shows  a  zone  of  losses  in  the  first  half  of  the  "supersonic 
diffuser"  characterized  by  a  great  pressure  increase  and  a  rapid  change 
into  a  subsonic  flow  with  an  adjacent  zone  of  increasing  velocity.  In  the 
narrowest  cross  section  of  the  diffuser,  velocity  of  sound  is  attained  the 
second  time.  Then  a  second  expansion  to  supersonic  velocity  takes  place." 
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The  detailed  discussion  of  his  experimental  results  is  condensed 
by  Ramm  in  the  following  theorems: 

"Theorem  1:  In  most  oases  the  diffuser  does  not  aot  as  a  devioe 
for  pressure  regain,  but  as  a  second  Laval  nozzle. 

"Theorem  2:  The  ratio  of  the  two  narrowest  oross  seotions  of  the 
diffuser  and  Laval  nozzle  determines  the  losses  in  the  diffuser. 

"Theorem  3*  The  transition  from  supersonio  flow  to  subsonio  flow 
in  the  "supersonio  diffuser"  does  not  ocour  by  a  simple  normal  compression 
shook. 


"Theorem  5<  1°  the  supersonio  diffuser,  pressure  losses  take  plaoe  which 

are  greater  than  those  of  the  normal  shock  wave  at  the  Mach  number  of  the 
nozzle.  The  losses  in  the  adjaoent  subsonio  diffuser  are  comparatively  small." 

Evidently,  this  summary  of  experimental  results  agrees  very  well 
with  the  important  qualitative  conclusions  of  the  present  theory. 

Due  to  the  restrictions  concerning  the  aoouraey  of  Ramses  tests  a 
numerical  comparison  of  the  data  lm  his  Fig.  10  is  not  given  in  general. 

Only  the  diffuser  efficiency  calculated  from  the  p*’  values  at  pressure 
equilibrium  and  at  Mach  Mo.  1.22  and  1.56  are  plotted  in  Fig.  20  (double 
circles) 1  those  at  Mach  Ho.  1.56  and  2.50  are  omitted  as  muoh  too  low. 

The  values  of  the  diffuser  efficiency  calculated  from  the  pressure  ratio 

of  the  entire  tunnel  (see  small  table  in  his  Fig.  10)  show  fair  agree* 
ment  with  the  more  exaot  measurements  at  the  same  tunnel  with  mercury 
manometers  (Fig.  20,  Curve  7),  but  are  omitted  due  to  larger  scattering. 

Values  of  F*»  at  pressure  equilibrium  from  Fig.  10(Ref.  3)  are  not  plotted 
in  Fig.  19.  They  agree  at  Mach  Ho.  1.22  with  the  other  measurements  at 
the  same  tunnel,  but  are  between  7%  and  15£  larger  at  the  other  3  Maoh  Nos. 

C.  Diffuser  Throat  Dimension  Measurements 


A  comparison  between  measured  and  calculated  diffuser  throat 
dimensions  (Fig.  19)  gives  a  very  reliable  oheok  of  the  theory,  due  to  the 
following  reasons: 

1.  The  value  of  the  diffuser  throat  depends  inversely  upon  the  value 
of  the  diffuser  efficiency  without  any  other  experimental  quantities 
involved.  (Equ.  36) • 

2.  The  value  of  F* •  according  to  (36)  is  connected  with  the  sonie 
velocity  in  P**  and  therefore  is  a  direct  oonsequenoe  of  the  primary  con¬ 
clusion  of  the  theory  that  subsonio  velocities  (^-solution)  takes  plaoe  in 
the  beginning  of  the  diffuser. 

3*  The  value  of  F**  is  Independent  of  all  losses  whioh  ocour  don- 
stream,  as  in  the  divergent  diffuser.  However,  those  losses  are  inoluded 
in  values  determined  by  the  total  pressure  ratio  of  the  entire  tunnel. 

(See  par.  D,  and  Fig.  20). 
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i;.  The  experimental  determination  of  F*1  is  simple.  It  represents 
automatically  an  average  across  the  area  and  does  not  involve  complicated 
pressure  measurements  which  are  very  sensitive  to  deviation  of  strictly 
one-dimensional  flow. 

5.  The  only  disadvantage  in  using  F**  is  a  possible  falsification 
of  its  value  by  the  boundary  layer  displacement  thickness.  Its  value  can 
be  scarcely  calculated  because  the  initial  thickness  at  the  end  of  the 
transformation  zone  is  unknown.  Nevertheless,  accelerated  subsonic  flow 
in  the  convergent  part  of  the  diffuser  will  keep  it  very  thin  as  it  is 
known  from  the  converging  part  of  the  Laval  nozzle. 

Fig.  19  shows  the  theoretical  curve  of  the  diffuser  throat  for 
f  *  0c6k,  according  to  Equ.  (36),  which  is  identical  with  the  corresponding 
curve  in  Fig.  9.  In  addition,  the  experimental  data  of  three  different 
series  are  plotted  from  tests  taken  in  the  supersonic  wind  tunnels  at 
Peenemuende  and  later  at  Kochel  (Ref.  16).  Test  Section  No.  1  and  2  had 
tunnels  of  the  same  size  and  the  same  diffuser  area  ratio  f  *  O.6I4.  but 
belonged  to  different  tunnels.  The  data  not  yet  published  are  taken  from 
test  records  and  graphs  in  the  hands  of  my  former  associate,  H.  U.  Eckert, 
who  carried  out  the  last  series  in  Feb.  19b5«  The  subsonic  values  are  taken 
from  an  average  curve  drawn  through  many  test  points.  The  subsonic  Mach  Nos. 
have  been  calculated  from  static  pressure  measurements  in  the  axis  of  the 
nozzle  by  means  of  a  long  probe.  The  F**  values  at  Ma  ■  J4..28  are  those 
obtained  after  supersonic  flow  through  the  tunnel  had  been  established. 

In  order  to  start  the  tunnel  approximately  10  -  20$  larger  F*’  values  were 
needed  producing  test  chamber  pressure  smaller  than  the  nozzle  exit  pressure. 

It  might  be  mentioned  that  the  second  throat  area  is  sensitive  to  inleakage. 
Sealing  the  tunnel  was  a  difficult  problem  for  Mach  number  3  and  higher. 

The  test  series  above  were  carried  out  with  glass  windows.  Other  tests  in 
February  19^5  with  the  so-called  steel  windows,  which  had,  as  it  was  known, 
insufficient  sealing,  needed  much  higher  F#*  values  (up  to  50$)  due  to  entrance 
of  additional  mass  and  increased  transformation  losses 0 

As  a  result.  Fig.  19  represents  an  excellent  quantitative  confirmation 
of  the  theoretical  calculations.  The  maximum  deviations  of  the  averaged  test 
points  from  the  theoretical  curve  are  +k%  and  -7$. 

D.  Experimental  Data  of  Total  Pressure  Ratio  of  the  Entire  Wind  Tunnel. 

There  are  test  data  available  to  the  author  concerning  the  total 
pressure  ratio  of  several  European  supersonic  wind  tunnels.  Assuming  that 
the  additional  losses  downstream  of  the  diffuser  throat  are  comparatively 
small  as  determined  experimentally  in  the  Peenemuende  tunnel  (see  Sec.  6,B, 
Theorem  5),  these  data  allow  a  third  check  of  the  theory,  completely  in¬ 
dependent  of  those  made  in  Par.  B  and  C.  of  this  section. 
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1«  r>inwiitnit  »  ltt»r  Koohel  Wind  Tunnel^ 


Th*  aaln  proof  In  Fig.  20  in  established  by  comparison  of  curve  A 
with  curve  7.  Curve  A  represents  tho  theoretical^- curve  for  f  •  0.61* 
as  shown  boforo  in  Pig.  6.  Curve  7  ■how  test  data  on  tha  Pc  enenu  ende-to  ohel 
fraa  jet  tunnol  (Hof.  16)  with  40  cm  x  40  on  nostle  axit  araa  and  50  cm «  5J  « 
diffuaar  intaka  araa,  corresponding  to  f  •  0.64.  Tha  taata,  not  yet  publiAed, 
are  carried  out  in  tunnel  Ho.  2  in  Koohel  by  wy  fomar  aaaoeiate  H.O.  Eokart 
and  raoordad  in  loohel  Graph  Ho.  PL  666  fro*  Mar  5.  *945 •  following 

oritloal  ranarka  ahould  be  Bade.  Purpoae  of  thaaa  taata  waa  tha  furnishing  or 
deaim  data,  and  not  ohaeking  of  any  theoretical  calculations*  Tha  central  aro 
shaped  nodal  support  waa  installed.  Due  to  tha  intermittent  tunnel  operation, 
tha  pressure  readings  (simple  mercury  manometers)  for  maximum  efficiency  were 
taken  just  before  breakdown  of  the  constant  conditions  in  the  teat  chamber,  the 
time  of  #iioh  sometimes  is  not  exactly  defined.  For  Mach  Ho.  4*3  the  tunnel 
needed  larger  throat  area  for  starting,  producing  pressure  in  the  teat  chamber 
smaller  than  the  no sale  exit  pressure  with  higher  losses.  This  throat  area 
oould  be  reduoed  during  operation.  It  is  uncertain  if  any  final  constant  state 
could  be  reaohed  during  the  short  running  time  of  about  16  seconds.  The  - 
value  plotted  in  ourve  7  corresponds  to  this  reduoed  area,  while  the 
belonging  to  the  first  area  would  lie  about  11*  below.  Hie  ^j-xalues  *?*••“* 
citive  to  the  exaot  test  section  Maoh  Ho.  for  Aich  the  KN8  i»  taken,  therefore 

average  value  in  the  tunnel  axis  over  at  least  ±  20  cm  was  used  for  evaluation. 
Considering  these  remarks,  the  agreement  between  theory  and  experiment  is  good, 
the  maximum  deviation  of  the  tests  is  about  6*  below  the  oalculated  ourve 
between  Ma  -  1  and  1.6.  This  oan  be  explained  by  additional  losses  in  the 
subsonic  diffuser.  The  snmll  positive  deviation  (2*)  at  Maoh  Ho.  4.3  i«  attri¬ 
buted  to  the  aforementioned  difficulties. 


Fig.  21  and  22  give  other  oowparisons  between  the  same  test  results  and 
the  theory,  using  other  variables  and  plotted  up  to  Mach  Ho.  10  in  a  semi- 
logarithmic  soale.  Fig.  21  shows  the  total  pressure  ratio,  the  test  points 
tltm  for  the  entire  tunnel  .  while  the  ourve  is  calculated  from  the 

dlffhser  theory-^.  Fig.  22  is  the  reciprocal  plot  of  Fig.  21,  showing  the  com¬ 
pression  ratio  needed  for  the  blowers  if  no  additional  losses  oocur.  This  repre¬ 
sentation  is  sometimes  favored,  for  instance  in  the  book  by  Liepmann  A  Puckett 
(Ref.  19)  Fig.  5.16.  The  curve  of  this  figure,  going  up  to  Mach  Ho.  4*2  and 
desi gnated  "conservative  values  based  on  performance  of  several  actual  wind 
tunnels,"  ie  praetioally  identical  with  the  test  values  of  the  loohel  tunnel 
show  in  Fig.  22.  Fig.  21  and  22  show  the  great  advantage  for  engineering 
purposes  of  having  a  diffuser  theory  extending  up  to  high  **«h  number. 

It  night  be  mentioned  that  P.  R.  Owen  (Ref.  20)  in  his  oareful  end  detailed 
siihi~t  on  equipment  and  test  results  of  the  Koohel  wind  tunnel  gives  oom- 
pres  si  on  ratios  in  the  range  from  Mah  Mo.  0.4  to  4*36  (Part  I,  Table  8),  Fhi  oh 

are  too  high  throughout  in  the  amount  of  5  to  20*.  Evidently  he  obtained  nie- 
1 coding  information. 
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It  was  the  author's  desire  to  prove  the  theoretical  results  in  other  free 
jet  wind  tunnel  data  than  those  from  Peenemuende  and  Kochel  but  this  is  possible 
only  with  limitation.  Data  are  available  from  the  Goettingen  AVA  and  the 
Aachen  AIA  free  jet  wind  tunnel 0 

2.  Goettingen  AVA  Wind  Tunnel 

The  Goettingen  AVA  free  jet  intermittent  type  supersonic  tunnel  has 
a  nozzle  exit  area  11  cm  x  13  cm  and  a  diffuser  intake  area  of  13  cm  x  15  cm; 
hence  an  area  ratio  of  0,75.  Tests  of  total  pressure  ratio  (see  Fig,  20, 

Curve  2)  were  carried  out  by  0.  Walchner,  but  are  reported  by  K.  Oswatitch 


The  test  results  lie  far  below  the  theoretical  curve  for  0.75,  which 
would  lie  somewhat  above  the  curve  A,  if  shown.  However,  the  diffuser  in 
this  tunnel  (Eef.  13)  has  a  completely  different  shape  from  that  generally 
used,  and  as  assumed  in  this  analysis  (Fig.  1).  it  has  only  the  convergent 
part  of  the  conventional  diffuser  up  to  the  minimum  area  and  then  expands 
abruptly  to  about  twice  that  area.  Thus,  a  kind  of  second  free  jet  with 
sonic  velocity  emerges  in  a  large  pipe.  Trying  to  calculate  this  process  by 
a  second  application  of  the  present  formulas  fails  and  does  not  explain  the 
low  diffuser  efficiency.  0.  Walchner  himself  designates  these  data  as 
preliminary  measurements  in  an  incomplete  war-time  diffuser"  (personal 
information  to  the  author).  The  data  themselves  cannot  be  checked  today. 

Dry  air  was  used.  K.  Oswatitsch  noticed  this  low  pressure  recovery  also, 
but  explained  it  as  belonging  to  free  jet  wind  tunnels  in  general,  which 
obviously  is  an  error. 

5o  Aachen  AIA  Wind  Tunnel  (10  cm  x  10  cm) 

The  Aachen  AIA  free  jet  intermittent  type  supersonic  tunnel  has  an 
exit  area  of  10  cm  x  10  cm  and  a  diffuser  intake  area  of  11  cm  x  13  cm,  hence 
an  area  ratio  of  f  »  0.70.  Under  the  direction  of  G.  Wieselsberger  (Ref.  15) 
the  author  built  and  tested  this  tunnel  from  193k  to  1937.  Test  records  (from 
Nov.  29,  1936)  concerning  the  total  pressure  ratio  of  the  entire  wind  tunnel 
with  model  installed  are  in  the  possession  of  the  author.  The  results,  not  yet 
published,  are  shown  in  Fig.  20  (Curve  6)  and  lie  considerably  below  the  theo¬ 
retical  curve  for  f  »  0.70  (not  shown).  The  difference  can  be  explained  by  the 
fact  that  humid  air  was  used.  Hence,  condensation  shocks  connected  with  total 
pressure  loss  and  change  in  Mach  number  and  velocity  distribution  will  have 
occurred,  all  of  which  are  functions  of  relative  and  absolute  humidity,  not  yet 
quite  known.  The  humidity  conditions  of  these  tests  are  likewise  unknown. 

Recent  evaluation  of  incomplete  tests  of  that  time  by  the  author  (Ref.  21)  show 
that  at  so-called  "medium  humidity  conditions"  (approximately  1+0  to  50#  relative 
humidity  at  Aachen)  condensation  shocks  occurred  with  a  total  pressure  ratio 
of  approximately  0.90  at  Ma  »  lJ+2  and  0.20  at  Ma  ■  2.22.  The  value  of  Ma  ■  3.15 
could  not  be  determined.  Using  these  values  for  evaluation,  the  test  data  lie 
then  in  fair  agreement  and  are  only  between  10%  and  k%  below  the  theoretical 
curve.  However,  this  agreement  cannot  be  considered  as  a  reliable  confirmation 
of  the  diffuser  thecry,  due  to  the  above  mentioned  uncertainties. 
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1*,  Aachen  AIA  Wind  Tunnel  (!<.  cm  x  12  cm) 


Investigations  of  the  flow  process  and  pressure  measurements  in  the 
Laval  nozzle  of  one  smaller  test  section  were  performed  by  H.  Eggink  (Ref 0  22) 
some  years  later*  One  part  of  his  tests  refers  to  a  closed  test  section, 
another  one  to  an  "open  jet  arrangement Whether  this  was  a  completely  open 
jet  or  a  "half  open  -  half  closed"  test  section  could  not  be  found  out  from  his 
report.  The  author  states  (page  20  of  the  translation)  that  a  "simple  estimate 
of  the  pressure  recovery  for  a  free  jet  is  not  possible  since  the  losses  in  the 
free  jet  are  not  yet  understood."  Nevertheless,  he  gives  in  his  graphs  37,  3*5 
a  comparison  of  the  pressure  recovery  in  the  closed  and  free  jet  test  section* 

The  term  "pressure  recovery"  used  in  title,  text,  and  figure  captions  referring 
to  closed  and  free  jet  test  section  is  misleading*  It  is  not  the  total  pressure 
ratio  as  comnonly  used  but  the  ratio  of  two  specific  static  pressures.  These 
pressures  have  not  been  measured  in  or  behind  the  diffuser  but  are  calculated 
out  of  pressure  measurements  in  the  Laval  nozzle  with  several  doubtful  assump¬ 
tions  for  one  specific  diffuser  expansion  ratio  behind  the  second  throat.  The 
report  does  not  contain  any  measured  static  or  total  pressure  in  or  behind  the 
diffuser.  Hence,  Eggink's  data  of  so-called  "pressure  recoveiy  of  the  open  jet" 
are  not  suitable  to  verify  the  present  diffuser  theory* 

Further  on,  Eggink  discusses  the  flow  structure  of  the  free  jet  arrange¬ 
ment  at  different  diffuser  openings  and  downstream  pressures  by  means  of  Schlier- 
en  pictures  in  his  figures  26,  27,  and  2g.  It  should  be  pointed  out  that  the 
diffuser,  which  he  comnonly  uses  is  extremely  short  and  steep  and  quite  different 
in  shape  from  the  slender  diffuser  type  (Fig,  1)  on  which  the  present  analytical 
treatment  is  based.  The  bell  mouth  of  his  diffuser  acts  at  the  same  time  as  the 
steeply  convergent  part*  Recapturing  cross  section  F2  and  second  throat  area 
F#’  are  identical.  No  space  for  a  transformation  zone  exists*  The  total  length 
from  diffuser  intake  lips  to  the  second  throat  is  0*625  H  only*  Therefore, 
quantative  agreement  between  his  experimental  results  and  the  present  theoretical 
calculations  cannot  be  expected.  But  it  is  of  interest  that  for  a  somewhat 
flatter  diffuser  intake  at  Ma  ■  1.3  and  pressure  equilibrium  in  the  test  chamber 
he  observes  (Fig,  26  g,  h  and  2&  g,  h)  a  strong  pressure  increase  and  transition 
from  supersonic  to  subsonic  speed  in  the  beginning  of  the  convergent  diffuser 
part.  This  fact  is  in  agreement  with  one  of  our  main  theoretical  conclusions* 
Eggink  interprets  it  falsely  as  a  normal  shock.  At  Ma  ■  2.0  no  actual  Schlieren 
pictures  are  given.  Schematic  sketches  showing  the  steep  diffuser  intake  and 
only  the  case  of.  test  chamber  pressure  smaller  than  the  nozzle  exit  pressure 
give  oblique  shocks  without  transition  to  subsonic  speeds. 

Eggink  himself  notices  (page  16)  that  the  flatly  converging  diffuser 
(Fig.  26  g,  h)  has  lower  losses  than  the  steeply  converging  diffuser.  (Fig, 

26  c,  d),  His  proposal  for  improvement  of  pressure  recovery  for  free  jet 
(Pg.  21,  Fig,  39)  ends  up  with  a  very  long  and  slender  diffuser  arrangement 
(distance  from  diffuser  lip  to  second  throat  is  approximately  5  H),  which  is 
similar  to  that  of  our  Fig.  1.  Nevertheless,  our  reasoning  is  quite  different 
from  his. 
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E.  Supplement:  Some  Wind  Tunnels  with  Closed  Test  Section 

Although  tunnels  with  closed  test  section  are  not  considered  in  this 
paper,  some  brief  remarks  in  connection  with  them  are  of  interest  for 
comparison.  For  this  purpose  *ig.  20  shows  also  the  diffuser  efficiency- 
calculated  from  total  pressure  data  of  some  supersonic  wind  tunnels  with 
closed  test  sections. 

The  data  of  the  Zuerich-ETH  wind  tunnel  with  continuous  operation  are 
given  by  J.  Ackeret  (Ref.  1+,  page  521)  and  involve  the  diffuser  losses 
only  without  losses  by  corners,  straightener  or  cooler  (Curve  No.  1).  The 
test  section  area  is  ItO  cm  x  38  cm.  The  pressure  distribution  measurements 
range  from  Ma  «*  1.00  to  Ma  =  2.75  and  have  been  extrapolated  by  Ackeret  up  to 
3.00.  The  air  was  not  dried  but  evidently  the  influence  of  humidity  is  low 
in  his  tests,  probably  due  to  reduced  relative  humidity  by  heating  up  the 
tunnel. 

The  data  of  the  Goettingen  KWI  intemittent  type  wind  tunnel  are  reported 
by  K.  Oswatitsch  (Ref.  lit)  in  a  special  study  on  pressure  recovery  in  the 
diffuser.  The  tunnel  dimensions  are  small,  6  cm  x  S  cm.  The  tunnel  was  used 
with  completely  closed  test  section  with  either  a  short  (31  cm)  cylindrical 
diffuser  (Curve  3)  or  a  long  (131  cm)  cylindrical  diffuser  (curve  1+)  and  some¬ 
times  with  half -open  -  half-closed  test  section  (curve  5)«  The  air  was  not 
dried  so  that  these  results  are  somewhat  affected  by  condensation  as  noted  by 
Oswatitsch  himself. 

Because  of  many  uncertainties  test  data  of  the  Guidonia  stipersonic  wind 
tunnel  with  closed  test  section  (1+0  cm  x  1+0  cm)  and  continuous  operation  are 
not  included  in  Fig.  20.  Gasperi  himself  (Ref.  23)  notices  that  in  the  Laval 
nozzle  ("Guidonia  Nozzle  No*  g"),  designed  for  adiabatic  expansion  to  Mach 
No.  2.11+,  large  losses  occur  which  lead  to  much  smaller  Mach  numbers 0  His 
measurements  of  total  pressure  inside  Idle  flow,  and  static  pressure  at  the 
tunnel  side  wall  do  not  fulfill  isentropic  relation  with  the  stagnation  pressure 
in  the  stilling  chamber.  Furthermore,  Gasperi  interprets  the  pressure  measured 
with  the  side  orifices  of  the  pitot  tube  as  static  pressure.  The  reliability 
of  this  interpretation  seems  doubtful.  Calculations  of  Mach  number  using 
different  combinations  of  these  pressures  from  his  Fig.  9  give  5  different 
values  of  Mach  numbers  ranging  between  2.07  and  1.6I+.  Gasperi  gave  the  values 
1,69  and  1.77  in  the  axis  of  the  flow  and  1.60  at  the  wall.  Similar  figures 
result  by  the  evaluation  of  his  Fig.  10.  This  fact  indicates  considerable 
losses  most  probably  caused  by  condensation  shocks  of  humid  air,  because  for 
dry  air  the  above  mentioned  isentropic  relation  between  the  pressures  has  been 
proved  to  be  correct.  Hence,  proper  evaluation  of  his  diffuser  pressure 
recovery  values  is  impossible.  It  is  unknown  to  the  author  if  later  measure¬ 
ments  exist.  In  any  case  the  value  Ma  “  2.00  given  by  Oswatitsch  (Ref.  11+) 
for  the  Mach  No.  of  this  nozzle  is  too  high. 

Not  included  either  in  Fig.  20  are  results  of  pressure  measurements  by 
W.  Kraus  (Ref.  2i+)  in  the  tunnel  No.  5  of  Peenemuende  with  16  cm  x  16  «cm  nozzle 
exit  area  and  half-closed — half -open  test  section.  In  19U+  test  chamber  and 
diffuser  were  brought  to  the  LFA  at  Braunschweig  and  connected  with  -the  com¬ 
pressor  system  of  the  supersonic  tunnel  there  for  continuous  operation.  Unfor¬ 
tunately,  these  tests  were  carried  out  with  humid  air  because  no  air  dryer  in  the 
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LFA  installation  was  available,  and  thus  all  results  are  influenced  consider- 
ably  by  condensation  shocks.  Efforts  have  been  made  by  the  author  to  evaluate 
these  tests  under  consideration  of  condensation,  but  without  success,  since 
some  data  are  missing  required  for  this  evaluation.  The  performance  of  these 
experiments  was  strongly  affected  by  war  time  events. 

Curve  8  represents  the  above  (Sec.  1,  C,  2)  mentioned  curve  established 
by  L.  Crocco  (Ref.  5,  Fig.  1)  and  designated  by  him  with  "pressure  recovery 
practically  obtainable  with  steady  compression."  It  is  converted  from  Crocco's 
definitions  r  into  the  S*  definition  used  in  this  paper  by  means  of  the 
corresponding  formulas  given  before  (Sec.  1,  C,  2  and  4)*  He  calls  the  curve 
a  "mean  experimental  value"  based  on  wind  tunnel  tests  by  Ackeret  between 
Ma  *  1  and  2.73  und  by  Castagna  at  Ma  “  4®2  (one  test  point).  he  trend  of 
Crocco's  curve  6  seems  to  be  fairly  optimistic  above  Ma  -  2.5  in  comparison 
with  the  closed  tunnel  data  of  Zuerich  ETH  (Curve  1)  and  Goettingen  KWI  (Curves 
3  and  4).  Its  fluctuations  have  no  real  meaning  but  are  caused  by  inaccurate 
readings  of  the  original  curve. 

F.  Comparison  Between  Diffuser  Efficiency  of  a  Free  Jet  and  a  Closed  Test 

Section  ^ind  Tunnel.  (Example) 

Fig.  20  shows  good  agreement  between  the  diffuser  efficiency  of  the  three 
closed  test  section  tunnels.  (Curves  1,  3,  4).  They  have  somewhat  higher 
diffuser  efficiency  (average  of  between  Mach  No.  1  to  3)  than  the  free  jet 
(f  “  0.64)  tunnel  (Curve  7).  The  efficiency  value  for  the  half-closed— half¬ 
open  test  section  tunnel  (Curve  5)  lies  somewhat  below  (average  of  11^)  that  of 
the  completely  free  jet  tunnel.  However,  it  would  be  incorrect  to  generalize 
these  conclusions.  For  a  correct  comparison  of  different  test  section  types 
one  has  to  consider  the  same  test  section  length.  As  shown  in  the  present 
analysis  in  case  of  the  free  jet  test  section,  the  diffuser  eff iciency  depends 

greatly  upon  the  area  ratio  f  (Fig.  6,  7,  2)  and  therefore  on  the  free  jet  test 
section  length  1^,  which  requires  this  area  ratio  f  (^ig.  IS).  The  energy 
ratio  decreases  with  increasing  test  section  length  (Fig,  17 ),  In  case  of  the 
closed  test  section,  no  such  analysis  is  known  up  to  now,  but  without  doubt  the 
diffuser  efficiency  will  likewise  decrease  with  Increasing  test  section  length. 

For  the  purpose  of  comparison  we  take  as  an  example  the  Zuerich  and  the 
Peenemuende-Kochel  tunnel,  both  of  practically  the  same  size  and  operating  at 
nearly  the  same  Reynolds  number  and  refer  their  efficiency  values  to  the  same 
test  section  length.  We  have  to  choose  the  test  section  length  of  l/H  ■  1.3 
because  this  length  was  actually  used  in  establishing  the  experimental  data 
of  the  closed  test  section  tunnel  at  Zuerich  (H  “  32  cm;  Curve  1  in  Fig.  20). 

The  pressure  recovery  data  of  the  free  jet  test  section  tunnel  at  Peenemuende- 
Kochel  were  measured  with  a  test  section  length  l/H  “  1.75  f  “  0.64.  In 
order  to  convert  these  experimental  data  to  l/H  “  1.3  we  find  first  liie  corres¬ 
ponding  area  ratio  f  “  0.74  from  Fig.  10  (f or  6 “  0.06l).  Then  using  the  present 
analysis  we  find  the  theoretical^,  belonging  to  f  “  0.74  from  Fig.  8  as  a 
function  of  Mach  number.  Finally  we  add  or  subtract  the  same  percentage 
deviation  from  this  theoretical^,  as  found  between  the  experimental  and  theo¬ 
retical  curve  for  f  ■  0.64  (difference  between  curve  7  and  curve  A  in  Fig.  20). 
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The  result  is  given  in  the  following  table: 


Comparison  between  diffuser  effioienoy  of  a  free  jet 
and  a  closed  test  seotion  tunnel  at  the  same  test 
section  length  l/R  »  1.3 


Converted  Experimental  Data 


Mach 

No. 

^FVee  Jet 

conv.  Koohel 
wind  tunnel 

^Closed  test  seotion  ^ 

Zuerioh  wind  tunnel  4 

^t>  Free  Jet 

^closed  Test  Section 

1.0 

0.859 

0.827 

1.039 

1.2 

0.810 

0.810 

1.000 

1.5 

0.779 

0.739 

0.937 

2.0 

0.768 

0.795 

0.966 

2.5 

0.763 

0.777 

0.982 

3.0 

0.753 

O.763 

0.995 

Remark: 

The  third  decimal 

Average 

has  limited  physioal  meaning. 

0.995 

It  is  seen  that  the  differenoe  in  diffuser  efficiency  between  both 
test  seotion  types  at  this  specific  length  is  very  small,  maximum  ±  at 
different  Mach  numbers  and  within  0.5#  as  average  in  the  given  supersonic  Uadi 
number  range  from  1  to  3,  Evidently  the  superiority  in  efficiency  of  the  closed 
test  seotion,  as  it  is  lcnoan  from  the  subsonic  range  (9.3.  see  Ref.  12)  is 
vanishing  in  the  supersonio  range,  according  to  this  example.  Between  Mach  So. 

1  and  2  there  might  be  still  a  slight  difference.  At  about  Mach  No.  3  both 
tunnel  types  have  praotioally  equal  diffuser  effioienoy. 

It  might  be  emphasised  that  even  though  the  decelerating  process  in  the 
diffuser  behind  the  olosed  test  section  will  be  quite  different  from  that  behind 
the  free  jet,  the  diffuser  efficiencies  of  both  test  section  types  are  nearly 
equal  in  the  supersonio  range  and  at  the  test  seotion  length  under  consideration. 
A  statement  for  Mach  Nos.  above  3  or  for  other  test  section  length  oan  be  made 
only  after  a  similar  theoretical  analysis  is  known  for  the  closed  test  section 
diffuser  as  it  is  given  here  for  the  free  jet  test  section  diffuser  or  after  more 
extended  experimental  data  are  established. 
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ITOi'ETCLATUES 


o.  »  Velocity  of  sound 

cP  «  Speoifio  heat  at  constant  pressure 

f  »  Ratio  of  nozzle  exit  cross  section  area  to  diffuser  intake  area 
F  ■  Cross  section  area  (if  with  subscript) 

F(x)*  Generalised  iaentropio  line 
g(x)»  Negative  first  derivative  of  FCx) 
g  •  Inverse  funotion  of  g 
6  »  Function  symbol 

H  ■  Height  of  nozzle  exit  cross  seotion 
K  *  Total  pressure  ratio 
l  »  Length  of  free  jet 
L  •  Energy 

m  m  Slope  of  isentropio  line 
M  ■  •  Dimensionless  velocity 

Ma»  •  Mach  Number 

p  *  Pressure 

po •  Total  pressure  after  losses 
P 

r  »  -po •  Dimensionless  pressure 
R  •  Gas  constant 
S  •  Entropy 

T  »  Absolute  temperature 
V m  Volume  flow 

W *  Compression  work  per  unit  mass 
w  ■  Velocity 

X  »  Velooity- variable  for  generalised  isentropio  lines 
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y  ■  Pressure* variable  for  generalised  i ■•ntropi o  lines 
(x,;y, )  «  Arbitrarily  oho sen  fixed  point 
(xrjVr)*  Point  of  oontaot 

y  •  Ratio  of  speoifio  heats;  y  ■  1.1*05  for  numerioal  calculations  of  air 
€  •  Half  of  expansion  angle  from  nozzle  exit  to  diffuser  intake  (see  Fig. 
C  •  Energy  ratio 
^  •  Efficiency  factor 
^  *  “f^s  (*••  •quation  (ll*)  ) 

0  •  Dimensionless  mass  flow 


9  •  Density 
«  Function  symbol 


SUBSCRIPTS 

*  Adiabatlo  prooess 
B  •  Blower 
C ■  Test  ohamber 
Cp.«  Abbreviation  for  Crooco 
D  •  Diffuser 
E  -  Energy 

Exit  of  diffuser!  intake  of  blower 
i  •  Inlet  of  diffuser 
IS*  Isotheralo  prooess 
kin»  Kinetio 
mox»  Maximum 
NS«  Kormal  shook 
P  •  Pressure 
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d •  Supersonic  solution 
(1  •  Subsonic  solution 
I  •  Houle  exit 

2-  Diffuser  intake  if  areat  properties  of  state  after  transformation 

process 

0«  Stagnation 

*  •  First  throat 

*  •  Second  throat 
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SCHEME  OF  SUPER  SONIC  WIND  TUNNEL 
WITH  FREE  JET  TESTSECTION 
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Fig.  / 

Scheme  of  supersonic  wind  tunnel  with  free  jet  test  section 
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